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ABSTRACT 

Aims. The abundance evolution of interstellar dust species originating from stellar sources and from condensation in molecular clouds 
in the local interstellar medium of the Milky Way is studied and the input of dust material to the Solar System is determined. 
Methods. A one-zone chemical evolution model of the Milky Way for the elemental composition of the disk combined with an 
evolution model for its interstellar dust component similar to that of Dwek (1998) is developed. The dust model considers dust-mass 
return from AGB stars as calculated from synthetic AGB models combined with models for dust condensation in stellar outflows. 
Supernova dust formation is included in a simple parameterized form which is gauged by observed abundances of presolar dust grains 
with supernova origin. For dust growth in the ISM a simple method is developed for coupling this with disk and dust evolution models. 
Results. The time evolution of the abundance of the following dust species is followed in the model: silicate, carbon, silicon carbide, 
and iron dust from AGB stars and from SNe as well as silicate, carbon, and iron dust grown in molecular clouds. It is shown that the 
interstellar dust population is dominated by dust accreted in molecular clouds; most of the dust material entering the Solar System at 
its formation does not show isotopic abundance anomalies of the refractory elements, i.e., inconspicuous isotopic abundances do not 
point to a Solar System origin of dust grains. The observed abundance ratios of presolar dust grains formed in SN ejecta and in AGB 
star outflows requires that for the ejecta from SNe the fraction of refractory elements condensed into dust is 0.15 for carbon dust and 
is quite small (~ 10~ 4 ) for other dust species. 
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1. Introduction 

In this paper we intend to study the population of dust grains in 
the interstellar medium of the Milky Way, its composition and 
evolution, and the input of such grains into the Solar System. 

Part of these grains are formed in stellar ejecta or stel- 
lar winds of highly evolved stars from the refractory elements 
therein and the resulting gas-dust-mixture is ultimately mixed 
with the general interstellar matter. But generally only some 
fraction of the refractory elements in the ejecta is really con- 
densed into solid phases; a big mass fraction of the refractory 
elements returned to the ISM - or sometimes even most of it - 
stays in the gas phase. Turbulent mixing in the ISM rapidly inter- 
mingles the material that is ejected by the many different stellar 
sources. While all the dust particles from stellar sources, called 
Stardust particles, retain their peculiar isotopic compositions of 
a number of elements indicative for their formation sites, in the 
gas phase the material from all sources is mixed together and 
the resulting isotopic composition of the mix is different from 
that of the Stardust particles. The refractory elements in the ISM 
gas therefore have different isotopic compositions than the same 
elements found in Stardust grains. 

Stardust grains are found in the Solar System as a rare 
fraction of the fine grained matrix material of meteorites (e.g. 
Bernatowicz & Zinner [19971 Hoppe 120041 Nguyen et al. p2TjfJ7l) . 
They are identified as such by the unusual isotopic composi- 
tion of at least one element which shows that the grains have 
condensed from material that contains freshly synthesized nu- 
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clei from stellar burning zones. Laboratory studies have found a 
big assortment of such Stardust grains, also called presolar dust 
particles, with a variety of chemical compositions that can be 
associated with a variety of stellar sources. The composition of 
the grain material indicates two basically different chemical en- 
vironments of formation, (1) a carbon rich environment, and (2) 
an oxygen rich environment, which yield two completely differ- 
ent groups of mineral compounds: 

1. solid graphitic carbon, diamond, silicon carbide, silicon ni- 
tride, and 

2. corundum, hibonite, spinel, magnesium-iron-silicates. 

Besides of these main components a number of minor compo- 
nents (e.g. titanium oxide, solid solutions of titanium carbide 
with zirconium and molybdenum carbide, kamacite, and cohen- 
ite) have been identified which are so far only known to exist as 
inclusions in grains of the major dust components. Additional 
components may exist and await identification, in particular 
since some kinds of dust may not survive all the stages between 
the stellar source and the final laboratory investigation. The ob- 
served isotopic anomalies indicate essentially two different kinds 
of stellar sources of the presolar dust grains, (1) AGB stars and 
(2) core collapse supernovae. 

From the instant of the formation in stellar ejecta on, star- 
dust grains are subject to destructive processes in the ISM by 
sputtering and shattering processes induced by supernova (SN) 
shocks (cf. Jones et al. 1996 and references therein). They are 
finally incorporated into newly formed stars and their planetary 
systems after about 2.5 Gyrs residence time in the ISM, which is 
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Table 1. Comparison of some observed properties of the galactic disk in the solar neighbourhood with model results 



Observable 


Model 


Observed 


Reference 


Total surface density X to t [Mqpc~ 2 ] 


56 


50-62 


Holmberg & Flvnnl2004l 


ISM surface density Ejsm [Mqpc~ 2 ] 


9.7 


7-13 


Dickevl 19931 






13 - 14 


Oiling & Merrifield|2001| 


Gas fraction Ei SM /S tot 


0.17 


0.05 - 0.2 




Surface density of visible stars X, [Mqpc~ 2 ] 


38.6 


30-40 


Gilmore et al.l 19891 


Surface density of stellar remnants [Mqpc~ 2 ] 


7.7 


2-4 


Meraetal.| 1998| 


Star formation rate B [MQpc~ 2 Gyr~'] 


3.1 


3.5-5 


Rana[T99T1 


SN II rate R S mi [pc~ 2 Gyr 1 ] 


0.016 


0.009 - 0.0326 


Tammann et al. 1994 


SN la rate /? SNIa [pc^Gyr 1 ] 


0.0024 


0.0015 - 0.0109 


Tammann et al.|l994| 


Infall rate [M pc~ 2 Gyr'] 


1.45 


0.5-5 


Braun & Thilker|2004l 



also the typical timescale for replenishment of the ISM with new 
Stardust. Theoretical studies have shown that typical lifetimes 
against destruction by SN shocks are of the order of only about 
0.5 Gyrs (Jones et al. 1996). This rather short timescale com- 
pared to the timescale for replenishment would result in a very 
low dust abundance in the ISM. This, however, is not observed. 
Instead one observes a high degree of depletion of the refractory 
elements in the gas phase of the ISM (e.g. Savage & Sembach 
1996; Jenkins [2004), and this clearly requires additional growth 
processes in the interstellar medium that tie up the atoms of the 
refractory elements in dust. The only possible sites where ac- 
cretion of gas phase material onto grains may proceed with a 
reasonable short timescale are the dense molecular clouds of the 
ISM (Draine 1990). Any solid phase material grown in molecu- 
lar clouds, the MC-grown dust, has isotopic compositions of the 
refractory elements different from that found in the Stardust par- 
ticles and may be identified by this property. Such dust material 
may be found both as coating of Stardust grains and as separate 
grains. 

Unfortunately, the isotopic composition of the refractory el- 
ements in MC-grown dust grains incorporated into the Solar 
System equals the isotopic composition of the elements in the 
Solar System, which makes it impossible to discriminate by lab- 
oratory investigations of isotopic abundances of refractory ele- 
ments alone between dust formed in the Solar System and MC- 
grown dust. There are other indications, however, which point 
to a presolar origin of some fraction of the interplanetary dust 
particles, the GEMS (Bradley [20031 Messenger et al. 12003]) , that 
shows isotopic abundances of refractory elements corresponding 
to normal Solar System isotopic abundances, and which there- 
fore are likely to be MC-grown dust grains. The category of 
presolar dust grains therefore includes also the MC-grown dust 
species which are isotopically inconspicuous, a property which 
makes them presently difficult to be identified as of extrasolar 
origin. 

If we intend to calculate theoretically the abundances of the 
different components of the interstellar dust mixture in the Milky 
Way at the solar cycle and in particular the composition of the 
dust mixture from which the Solar System formed, we have to 
construct a model of the Milky Way's chemical evolution that is 
coupled with a model for the evolution of the dust component of 
the interstellar matter. The evolution of the dust component is not 
independent of the evolution of the element abundances since the 
refractory elements forming the dust are only gradually formed 
during the course of the chemical evolution of the galaxy. The 
model for the dust evolution needs to consider the injection of 
Stardust into the ISM, the destruction processes of dust in the 
ISM, and the growth processes in molecular clouds. Very simple 



models for the evolution of the dust content of galaxies have al- 
ready been constructed (e.g. Lisenfeld & Ferrara 1998; Hirashita 
2000; Edmunds HED Morgan & Edmunds |2003j Inoue[2003]), 
but these are too simplistic to allow for a detailed calculation of 
the composition of the interstellar dust mixture. Only the method 
developed by Dwek ( 1998 ) to integrate the chemical evolution of 
the galactic disk and the dust evolution into a common model is 
sufficiently detailed to allow a modeling of the complex inter- 
play between the processes determining the dust evolution and 
has the potential of being extensible to even more complex sys- 
tems. This model is a one-zone model, i.e., the galactic disk is 
approximated by a set of independent cylinders with all physical 
variables within a cylinder averaged over the vertical direction 
with respect to the disks midplane, and a one-phase model, i.e., 
one averages the properties of the ISM over its different phases 
(cold, warm, and hot; cf. Tielens 2005). This type of model al- 
lows a successful and at the same time rather easy calculation of 
some important properties of the Milky Way disk, in particular 
of its chemical evolution (cf. Matteuchi 2003 ). The price one has 
to pay for the simplifications is that some processes, in particular 
those depending critically on the phase structure of the ISM, can- 
not be treated with sufficient accuracy. Nevertheless, the results 
obtained by Dwek (1998]) show that such a simple model can 
successfully be used to calculate the evolution of the interstellar 
dust. We take the model of Dwek as a basis for constructing a 
model which allows to treat a more complex mixture of Stardust 
and MC-grown dust. 

For the input of Stardust from AGB-stars we use our re- 
cent results for the dust production by AGB-stars (Ferrarotti & 
Gail|2006 ), which are somewhat extended. These tables present 
rather detailed information on the amount and composition of 
Stardust formed by AGB-stars. For Stardust from SN we follow 
the procedure of Dwek ( 1998 ) and use a simple parameterization 
for the dust production rate since no suited other information on 
dust production is available. Observations are inconclusive and 
theory is only in its infancy (Schneider et al. 2004; Nozawa et 
al. 120031 and references therein). We, in turn, try to gain some 
insight into the dust production efficiency of supernovae by com- 
paring our model results with meteoritic abundances of Stardust 
from SNe. The dust growth in molecular clouds is treated in 
more detail as in the model of Dwek ( 1998 ) since we intend to 
discriminate in the model between Stardust and MC-grown dust. 
In principle it would also be necessary to consider that according 
to observations of element depletion in the ISM the MC-grown 
dust has a definite core mantle structure with a more resilient 
core and a more easy destructible mantle. A theoretical treat- 
ment of this grain structure would require to consider at least a 
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two-phase interstellar medium (cf. Tielens 1998J Inoue[2003), 
and not a simple one-phase model as in our present calculation. 

The plan of our paper is as follows: In Sect. [2]our evolution 
model for the solar neighbourhood of the Milky Way is intro- 
duced and some results for the chemical evolution are discussed. 
In Sect. [3] the model for the dust return by stars is explained. 
Section|4]discusses the dust destruction and growth processes in 
the ISM and Sect.[5]presents the results for the evolution of the 
interstellar dust. Some concluding remarks are given in Sect. [6] 

2. Chemical evolution 

To study the evolution of the dust content of our Galaxy we de- 
velop a standard open model of galactic chemical evolution. In 
this model the Milky Way is formed by slow infall of primor- 
dial gas from the halo or intergalactic space. Merging with other 
galaxies seems not to have played a mayor role during most of 
the lifetime of the Milky Way, except for the very first evolu- 
tionary phase for which stellar dynamics (cf. Helmi et al[2006) 
and elemental abundances (Reddy et al. 2006; Ramirez et al. 
120071 1 indicate that there were major merging events. Merging 
is not considered in the model. In the one-zone approximation 
we neglect radial motions in the galactic disk (but cf. Vorobyov 
& Shchekinov 2006 ) and consider its evolution in a set of inde- 
pendent rings. 

We solve numerically a classical set of non-linear integro- 
differential equations for the chemical evolution of the Milky 
Way following a mathematical formulation similar to Dwek 
( 1998 ), who first extended the standard system of equations for 
the chemical evolution of the galactic disk (cf. Matteucci 2003 ) 
to include the evolution of the dust component of the Galaxy. 
The basic set of equations for the surface densities of gas, stars, 
nuclei, etc. is not repeated here but can be found in the papers 
cited before. We specify in the following only our choices for 
some important input quantities for the model calculations. 

2.1. Basic model parameters 
2.1.1. Infall 

First models of chemical evolution of the Galaxy were simple 
"closed-box" models, in which the Galaxy's mass is already 
fixed at the initial instant of evolution. However, these models 
failed to reproduce the metallicity distribution of metal poor 
stars, one of the most important observational constraints on 
chemical evolution modeling, a problem being known as G- 
dwarf problem. Open models assume formation of the Galaxy 
by accretion of primordial or metal poor gas from extragalac- 
tic sources to solve the G-dwarf problem, as was first suggested 
by Chiosi (1980), and later discussed by Pagel (1997). In open 
models the total surface density of the disk changes by the accre- 
tion of gas, outflows, and radial motions within the disk. For the 
Milky Way outflows can be neglected due to the strong gravita- 
tional potential, and radial motions are neglected in the one-zone 
approximation. In our model the infall rate entirely defines the 
evolution of the total surface density. 

Several scenarios for gas accretion have been proposed by 
different authors, suggesting different rates and sequences of for- 
mation of the galactic components, see Matteucci (2003) for de- 
tails. Models assuming an exponentially decreasing infall of the 
gas are most successful in reproducing the G-dwarf distribution. 
Following Chiappini et al. d 19971 ) we adopt a two-infall exponen- 
tially decreasing model that assumes two subsequent episodes of 
Galaxy formation. Initially, the halo and thick disk are formed 



during a short period of about th ~ 1 Gyr, then the thin disk is 
formed by accretion of material on a much longer timescale of 
td ~ 7 Gyr at Solar galactocentric radius (here r© = 8.5 kpc). 
The accretion rate is given in this model by the expression: 



dSi(r,f) 



inf 



df 

(XDinfAMe-^H 
(X,)infA(r)e- f ^ 



for t < f thin 

(XXMr^'-'^v for t > ftMn 



(1) 



where f t hm = 1 Gyr is the time of onset of accretion onto the thin 
disk. The formation of the disk is assumed to start ?q = 13 Gyrs 
ago. 

The coefficients A{r) and B(r) are derived such as to repro- 
duce the present day density distribution of the disk. At the solar 
cycle one has 



A(r Q ) 
B(r Q ) 



r H (l-e-'o/TH) 

-''tot 



r D (l 



- ('G -? lhin)/ T D 



(2) 
(3) 



For details we refer to Chiappini et al. (1997), Alibes et al. 
(2001 ). The value T. (ot (rQ, ?g) for the current total density of the 
disk is taken to be 56 Mq pc~ 2 according to Holmberg & Flynn 
(2004). For the contribution of the thick disk to the total surface 
density we choose "Eu = lOM0pc~ 2 . (Xi)i„f and (X;)-„y denote 
element abundances of infalling gas, which we assume to be pri- 
mordial. 



2.1.2. Stellar birthrate 

Observations of global star formation rate in spiral galaxies sug- 
gest a Schmidt-law type of dependence of the stellar birthrate 
on some power of the total gas surface density B oc E" as with 
n » 1.5 (Kennicutt 1998). An additional dependence of the star 
formation rate on the total surface density 2 tot (r, t) was suggested 
in self regulating star formation theory (Talbot & Arnett [19751 ). 
Later Dopita & Ryder (1994 ) confirmed this by observations and 
suggested an empirical law of star formation B(r, t) oc 2" ot E'" 
with m = 5/3 and n = 1/3 giving the best fit for the observed re- 
lationship between the stellar brightness and the surface bright- 
ness in H a in galactic disks. We adopt these values for the powers 
in the star formation law and choose analogous to Alibes et al. 
(120011 ) the following form of the star formation rate: 



S tot (r,Q"Z g (r,Q" 
B(r, t) = v 



xn+m—l 



(4) 



We take into account a star formation threshold, i.e. a minimum 
surface density S g required for star formation, which is set to 
7 M0pc~ 2 (Kennicutt 1998). In the numerical calculation the 
transition from zero to the threshold value is smoothed in order 
to avoid the unphysical numerical oscillations in the solution for 
the surface density S g close to the threshold that are produced 
otherwise by some integration methods. The constant v is fitted 
such that the model fits the present day star formation rate in the 
solar neighbourhood. We take v = 1.3 Gyr 1 . 

2.1.3. Stellar lifetimes 

We refuse the instantaneous recycling approximation, i.e. the as- 
sumption that massive stars die immediately after their birth and 
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return metals to the ISM, and consider stellar lifetimes as a func- 
tion of stellar mass and metallicity using an analytical approxi- 
mation given by Reiteri et al. d!9961 l. The formula of Reiteri et 
al. ([1996 ) is a good fit for the stellar lifetimes computed by the 
Padova group (Alongi et al. 1993; Bressan et al. 119931 Bertelli 
et al. 1 19941 ) in the metallicity range 7 • 10" 5 < Z < 3 • 1(T 2 and 
for initial masses between 0.6 and 120 Mq. 

2.1.4. Initial mass function 

The initial mass function (IMF) O(M) gives the distribution of 
stellar masses born in a star formation event. The IMF is one of 
the most important components of the galactic chemical evolu- 
tion model as it establishes the frequency of low and high mass 
stars, and thus their relative role in the chemical evolution. The 
IMF is commonly assumed to be constant in space and time and 
is usually approximated by a power law <I>(M) = AM~ a , where 
A is a constant derived from a normalization of the IMF to unity 
in the considered mass interval {M/,M„}. The first such an IMF 
for Solar neighbourhood stars was proposed by Salpeter ( 119551 
with a power law index a = 2.35, which is still widely used in 
chemical evolution models. Later studies showed that the IMF 
is signif icantly flatter below 0.5 M Q (e.g. Miller & Scalo |1979l 
Kroupa [19931 1. For a discussion of the concept of an IMF see 
Scalo (120051 and Elmegreen & Scalo ( |2006b . 

In the present paper we adopt a frequency distribution of stel- 
lar masses at star formation consisting of four separate power- 
law type distributions in four separate intervals of initial masses 
proposed by Kroupa (20021: 



<D(M) = A 



CiAT 03 , 0.01 < M/Mq < 0.08 
C 2 AT L3 , 0.08 < M/Mq < 0.5 
C 3 AT 2 - 3 , 0.5 < M/Mq < 1.0 
C 4 M~ 2 - 7 , 1.0 < M/Mq < 100, 



(5) 



where the coefficients C\ = 2.0158, C 2 = 0.1612, C 3 = C 4 = 
0.0806 are derived from the normalization procedure. Masses 
are in solar masses. During the calculation we found, that one 
obtains better model fits if for high mass stars the exponent is 
changed to 2.55. Such a somewhat flatter power law (cD oc M~ 26 ) 
is, for instance, observed for massive stars in the Orion nebula 
(Preibisch et al. 120021 

The average stellar mass is then given by the integration over 
the full range of stellar masses: 



M,, 



1 

0.01 



MO(Af) dM 



(6) 



2.1.5. Nucleosynthesis Prescriptions 

Nucleosynthesis prescriptions and stellar yields are another im- 
portant ingredient of chemical evolution modeling. Usually stel- 
lar yields are divided into three main categories according to 
stellar masses, and result from extensive stellar evolution cal- 
culations. 

The single low and intermediate mass stars from the mass 
range 0.8-8 M0 contribute to the enrichment of the Milky Way 
with heavy elements due to excessive mass-loss during the fi- 
nal stage of their AGB evolution. We adopt the yields for H, 
4 He, 12 C, 13 C, 14 N, and lfi O from van den Hoek & Groenewegen 
(1997), tabulated for the range 0.8 - 8 Mq of initial masses 
and for metallicities from 10~ 3 to 4 ■ 10~ 2 . For 23 Na, 24 Mg, 
25 Mg, 26 Mg, 26 Al, and 27 Al the yields of Karakas et al. ([20031 ) 
for the mass-range 1.0 - 6.5 Mq and range of metallicities 



Z = 0.004, 0.008, 0.02 are used. Outside of the range of tables 
the data are extrapolated. 

Rates of SN la explosions are calculated in the approxima- 
tion of Matteucci & Greggio ( 1986 ), based on the classical sce- 
nario of deflagration in C-0 White Dwarfs in binary systems 
(Whelan & Iben 119731 ), with a modification recently proposed 
by Hachisu et al. ( 119961 119991 ) that accounts for the important 
effect of metallicity on mass transfer in binaries to a compact 
object. In this scenario the accreting White Dwarf develops an 
optically thick wind. If [Fe/H] < -1 the wind is too weak for 
a SN la to occur; only for higher metallicities one has a contri- 
bution from SN la explosions to the heavy element production. 
The SN la yields are taken from Iwamoto et al. ( 1999 ). The pa- 
rameter determining the frequency of events is fitted such that 
the iron abundance of the Solar System is reproduced; a value of 
P — 2 x 10~ 2 was found to give the best result. 

The problem of yields from massive stars is complicated by 
the necessity to model supernova explosions, many details of 
which are still unknown. We adopt the recent nucleosynthesis 
prescriptions by Nomoto et al. ( 120061 ) which presented SN II nu- 
cleosynthesis yields as a functions of stellar masses (from the 
range of 1 1-50 Mq), metallicity, and explosion energy. Two dis- 
tinct new classes of massive SNe are taken into account in this 
work: very energetic Hypernovae, and very faint and low energy 
SNe. Nucleosynthesis in hypernovae can explain the observed 
trends of (Zn, Co, V, Ti, Mn, Cr)/Fe in extremely metal poor stars 
(Nomoto et al. [20061 Kobayashi et al. 120061 ). For comparison 
purposes we implemented also the yields of Woosley & Weaver 
J 1995b for stars with masses from the range of 11-40 Mq, that 
are most commonly used in chemical evolution calculations. 

For massive stars with M > 40 Mq the mass returned by the 
stars up to the end of carbon burning is taken from the models 
of Sch aller e t al. ( [19931 , Schaerer et al. (fT993l , and Charbonnel 
et al. d 1 9931 >. It is assumed that the remaining mass collapses 
into a Black Hole. The mass-return of nuclei is determined from 
the models for all those nuclei, for which surface abundances are 
given in the tables. For all others we assume that their abundance 
in the returned mass equals their initial abundance. 

It is assumed that the mass returned by stars is mixed with 
the general interstellar medium on timescales much shorter than 
the timescale for conversion into new stars, i.e., the composition 
of the interstellar medium is assumed to be homogeneous at each 
instant. This is justified by the observed low scatter of element 
abundances in the present ISM and of stellar element abundances 
in open stellar clusters (see Scalo & Elmegreen 2004, and refer- 
ences therein). 

2.2. Evolution of some disk properties 

We now show some results of a numerical calculation of the 
galactic disk's chemical evolution at the solar cycle that are im- 
portant for our problem. For the model presented in the follow- 
ing the nucleosynthetic yields of massive stars are taken from 
the tables of Nomoto et al. (120061 ). 

Figure Q] shows the evolution the total surface mass density 
£tot( r O' ^ anc ^ °^ t ^ le interstellar medium Hism^q, f) for the 
galactic disk. In the model it is assumed that the formation of the 
disk started 13 Gyrs ago. Initially most of the material in the disk 
was in gaseous interstellar matter, today and at the time of Solar 
System formation most of the galactic matter is condensed into 
stars. A minor fraction is locked up in stellar remnants (White 
Dwarfs, Neutron Stars, Black Holes). 

The evolution of the stellar birthrate is shown in Fig. [2^. Star 
formation commences about 1 Gyr after the onset of matter in- 
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Fig. 1. Time evolution of the total surface density E tot (full line), 
of the surface density of visible stars (dotted line), and of the in- 
terstellar matter surface density Eism (dashed line) of the galactic 
disk at the solar distance from the galactic centre, and observed 
values for the present day total surface density in the solar neigh- 
bourhood (Holmberg & Flynn 20Q4)), for the surface density of 
the stellar component (Gilmore et al. 1989), and for the surface 
density of the interstellar medium (Dickey [1993l l. 



fall since about 1 Gyr time is required in the two-infall model 
until the gas density at the galactocentric distance of the sun in- 
creases to the threshold value for star formation of Mism = 7 M0 
(Kennicutt |1998l l. The stellar birthrate culminated about 2 Gyrs 
after the onset of star formation and since then it gradually de- 
clines. Most of the stars born are low and intermediate mass 
stars; the massive stars mass fraction of the newly born stars is 
only 6.5% according to the initial mass function Eq. ©, but this 
small fraction is responsible for nearly all of the heavy nuclei 
synthesized in the Milky Way. 

Figure [2J? shows the evolution of the supernova rates at the 
solar galactocentric distance tq. Because of the short lifetime of 
massive stars the supernova rates for type II supernovae closely 
resembles the birthrate of stars. Supernovae of type la appear 
with a delay of several Gyrs because (i) their progenitors are 
long lived intermediate mass stars and (ii) the suppression of 
supernova explosions in binaries at low metallicities proposed by 
Hachisu et al. ( 19961 119991 . Since supernovae of type la are the 
main sources of Fe, the iron abundance increases in the Milky 
Way only on a rather long timescale. 



2.3. Chemical evolution of the disk 

The viability of galactic chemical evolution models is usually 
tested by comparison with some standard observational con- 
straints: G-dwarf metallicity distribution, age-metallicity rela- 
tion, Solar System abundances at the instant of its formation t$$p, 
and evolution of element abundance ratios over time. 




Fig. 2. (a) Evolution of the astration rate B at the solar galacto- 
centric distance tq. The errorbar shows the presently observed 
stellar birthrate (Rana fT99TT l. (b) Evolution of the supernova type 
II (full line) and type la (dashed line) rates at the solar galac- 
tocentric distance rQ and observed values at present time from 
Tammann ( f!994b . 



2.3.1. Evolution of metallicity 

Figure [3] shows the time evolution of the metallicity Z and the 
abundance ratio [Fe/Hj] of the interstellar medium at the solar 
radius rQ. The predicted evolution of the [Fe/H] ratio in the ISM 
is compared with the age-metallicity relation of the solar neigh- 
bourhood of late-type dwarfs from Rocha-Pinto (2000). The thin 
black vertical line shows the instant of Solar System formation 
?ssf, which we assume to be 4.56 Gyr ago. The filled circles indi- 
cate the observed Solar metallicity Zq and the present day ISM 
metallicity. The model reproduces the observational values quite 
well. 

All model calculations for the evolution of heavy element 
abundances with time predict a well defined relation between 
metallicity and time-of-birth at a certain location in the galactic 
disk like that shown in Fig. [3] Observationally determined ages 



The abundance ratio [X/Y] for two elements X and Y is defined as 



[X/Y]=log(g)-log(gj , 



O 

where ex is the element abundance of element X by number relative to 
hydrogen. 
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Fig. 3. Evolution of metallicity Z of the ISM at the solar galac- 
tocentric distance rQ and of the [Fe/H] abundance ratio. The er- 
ror bars show the observed age-metallicity relation from Rocha- 
Pinto (2000). The thin vertical line indicates the Solar System 
birth time, and the two filled circles indicate the observed metal- 
licity of the sun and of the present day ISM 




Fig. 4. G-dwarf metallicity distribution in the solar vicinity pre- 
dicted by the model and the observed distribution as derived 
by Nordstrom (2004). The thin dashed line shows the G-dwarf 
distribution from direct calculations, while the thick dotted line 
is the result of a convolution with a Gaussian with halfwidth 
0.2 dex to account for the observational scatter. 



obtained by comparing the position of a star in the Hertzsprung- 
Russel diagram with evolutionary isochrones, and relating spec- 
troscopically determined metallicities of stars with such age de- 
terminations, show a tremendous scattering of metallicities for a 
given age. It has been concluded that this reflects (i) a true scat- 
tering of metallicities of the matter out of which stars are formed 
at given galactocentric radius and birthtime, and (ii) possibly 
a mixing of stars from different galactic zones by radial diffu- 
sion (Edvardsson et al. 1993). Pont & Eyer (2004) have shown, 
however, that the tremendous scattering results most likely from 
the difficulty of obtaining reliable stellar ages from evolutionary 
isochrones and that any true internal scattering of metallicities at 
given age is probably less than 0.15 dex. More careful analysis 
of age-metallicity relations based on such improved methods (da 
Silva et al. 120061 ) also seem to support a small intrinsic scatter- 
ing of metallicities at a given age. The age-metallicity relation 
of Rocha-Pinto (2000 ) is based on ages determined from spec- 
troscopic indicators of chromospheric activity, which cannot be 
considered as particular accurate, but also this points to a small 
internal scattering of metallicities. 

The age-metallicity relation is reasonably well reproduced 
for the last about 10 Gyrs, but there is an increasing discrepancy 
for earlier times. This is a general problem of all such evolution 
calculations and results most likely (i) from the assumption that 
the infalling material has primordial abundances while in fact 
there is some pre-enrichment with heavy elements in this mate- 
rial, and (ii) from unrealistically high stellar ages for many stars 
due to the rather crude methods of age determination. 

2.3.2. Metallicity distribution of G dwarfs 

The G-dwarf metallicity distribution is one of the most important 
observational tests, since lifetimes of G-dwarf stars are compa- 
rable with the Galactic age, so that their metallicity distribution 
reflects past star formation history. We compare the observed 
G-dwarf distribution from the most recent and most complete 
compilation of Nordstrom et al. (2004) with that predicted by 
the model in Fig. [4] The dotted line shows the direct calcula- 
tional results, while the dashed line displays a convolution of the 



model results with a Gaussian with dispersion of 0.2 dex in order 
to simulate observational errors in the metallicity determination 
and intrinsic cosmic scatter in metal abundances. The errors of 
modern abundance determinations are usually claimed to be 0. 1 
dex or even less. The true scatter of stellar abundances for stars 
born at the same instant and location is difficult to determine 
since for single stars neither their birthplace nor their birthtime 
is accurately known. The small scatter of abundances between 
stars in open stellar clusters indicate, however, that the intrinsic 
scatter seems to be very small (see Scalo & Elmegreen 2004 
and references therein); we arbitrarily assume a contribution of 
0. 1 dex to the total scatter. 

The calculated metallicity distribution reproduces the gen- 
eral trends of the observed distribution, but it does not agree 
particular well. Significant deviations are seen for low and high 
metallicities. After convolution the discrepancies at the higher 
metallicity end disappear almost completely. The theoretical re- 
sult then reasonably well agrees with the observations. For low 
metallicities the discrepancies persist and indicate that our model 
assumptions are likely not realistic for the earliest evolutionary 
phase. Since for the main application of our model this phase is 
not important, we did not try to improve the model in this re- 
spect. 

2.3.3. Evolution of abundances of individual elements 

Abundance ratios [X/Fe] of elements X and their variation with 
time reflect the synthesis of heavy elements during galactic evo- 
lution. The reproduction of these variations by the model is one 
of the most important tests for the reliability of the model. For 
comparing the variation of [X/Fe] with observed variations of 
stellar abundances, stellar ages would be required which, how- 
ever, are unknown or of low accuracy. One prefers to compare in- 
stead the variation of the abundance ratios [X/Fe] with the abun- 
dance ratio [Fe/H], since [Fe/H] is also determined from stel- 
lar atmosphere analysis and varies, at least for the Milky Way, 
monotonously with age of the galactic disk (cf. Fig. [3), i.e., can 
be taken as a measure of stellar age. We have calculated in our 
model the evolution of 63 isotopes using nucleosynthesis pre- 
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Fig. 5. Comparison of the predicted abundance ratios of main dust forming elements [El/Fe] with observations of stellar abundances. 
The solid and dashed lines show model calculations with Nomoto (2006 ) and Woosley & Weaver ( 1995) SNII yields, respectively. 
We corrected WW95 yields for Fe and Mg to achieve better fits to observations. For illustrative purposes a model calculation with 
uncorrected Mg yields from WW95 is shown with a thin dashed line. The observed stellar element abundances for F and G stars 
from the solar neighbourhood are shown with different symbols for each of the sources (Akerman et al. 2004, Reddy et al. 2003 
Soubi ran et al.|20 5|] Me lendez et al . 120021 Jonsell et al. l2005l Venn et al. l2004l Chen et al. l2000l Gratton et al. fTWTl Caffau et al. 
|2005]Cayrel e t al. 120041 



scription of Nomoto et al. ((2006) and, for comparison, that of 
Woosley & Weaver d 1 995b . Results are presented in Fig. [5]for 
the elements related to dust formation. We concentrate here on 
these elements, since in the present work we are mainly con- 
cerned with problems related to interstellar dust evolution. 

The figure shows as thick lines the model results if SN II 
yields from Nomoto et al. ([2006 ) are used, and as thick dashed 
lines the corresponding results if yields from Woosley & Weaver 
(1995 ) are used. The various dots, crosses etc. show results of 
stellar abundance analysis for G stars from the solar neighbour- 
hood; the sources of data are given in the figure caption. These 
data show a considerable scatter because of the errors of abun- 
dance determinations and possibly some small intrinsic scatter 
of element abundances of stars of comparable age. Nevertheless 
there are clear observable correlations between the abundance 
ratios [X/Fe] and [Fe/H]. For the elements shown the new re- 
sults of Nomoto et al. give better agreement between the cal- 
culated abundance evolution and the observed correlations of 
[X/Fe] with [Fe/H] than the older Woosley & Weaver results, 
for other elements, however, there are some discrepancies with 
observations. 

With the yields of Woosley & Weaver ( 1995 ) there are some 
substantial problems. First, the iron yields of Woosley & Weaver 
are too high, as already found in Timmes et al. ( 1995), and we 
follow their recommendation to reduce the Fe yield. Second, 
there is another severe problem with the Woosley & Weaver re- 
sults for magnesium. The calculated abundances based on the 
original yields are shown in Fig. [5] with a thin dashed line. These 
abundances are definitely too low, a problem which is known 
since long (e.g. Goswami & Prantzos 2000; Francois et al. 2004). 



A comparison with the observed evolution of stellar magnesium 
abundances with metallicity shows that the shape of the [Mg/H]- 
[Fe/H] -relation is reasonably well reproduced by the model, ex- 
cept that the absolute values of [Mg/H] are systematically too 
low by a factor of 2.5. We have increased therefore the Mg yields 
of Woosley & Weaver ( 1 995 ) by this factor in order to reproduce 
the Mg abundance of the Solar System. The resulting variation 
of [Mg/H] with [Fe/H] is shown in the figure with a dashed line, 
which reproduces the observations much better. Such a correc- 
tion would be necessary for the purpose of calculating dust abun- 
dances, since reliable results for dust condensation require that 
the abundance ratios Si/Mg and Fe/Si of the main dust-forming 
elements agree with the observed abundance ratios in the Milky 
Way. Otherwise one would get a deviating dust mixture. 

Since the yields of Nomoto et al. (2006) gives results for 
the abundance evolution of the main dust forming refractory el- 
ements much closer to observations than the Woosley & Weaver 
( 1 19951 ) yields, and since they do not require for this to introduce 
some ad hoc scalings, we prefer to use the Nomoto et al. (2006) 
yields for the model calculations. 

Figure [6] compares the (adjusted) calculated abundance ra- 
tio Si/Mg with observed abundance ratios in the atmospheres 
of nearby F and G stars and their correlation with metallicity 
[Fe/H]. The model results are close to the observed values. The 
Si/Mg ratio determines the nature of the silicate dust that can 
be formed if Mg and Si both are completely condensed into 
dust. Oxygen is in any case abundant enough for formation of 
any kind of Mg-Si-compound. For a ratio of Si/Mg = 1 one can 
form enstatite (MgSiOs), for Si/Mg = 2 one can form forsterite 
(Mg2Si04). For values in between a mixture of both can be 
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Fig. 6. Abundance ratio Si/Mg of the major silicate dust forming 
elements. The full line corresponds to a model using SN yields 
of Nomoto et al. (2006 ), the dashed line to a model using SN 
yields from Woosley & Weaver (1995). For the latter the Mg 
abundance is scaled such that it reproduces the solar Mg abun- 
dance at [Fe/H] =0. 



formed; here part of the Mg can be replaced by Fe and one 
rather forms a mixture of magnesium-iron-silicates. The stellar 
data show that the Si/Mg ratio is close to unity for metallicities 
[Fe/H] > -1, and therefore one observes the formation of Mg- 
Fe-silicates in space. 

2.3.4. Solar System abundances 

The model should also reproduce the element abundances of the 
Solar System since they reflect the ISM composition at at the 
instant fssF = 4.56 Gyr ago, when the Sun was formed. Table [2] 
shows element abundances in the Solar System in the frequently 
used logarithmic scale (e is the abundance of an element relative 
to H by number) 



a = log e + 12 



(7) 



for the elements from H to Zn that can be compared to the re- 
sults of the model calculation if we use supernova yields from 
Woosley & Weaver ( 1 19951 ) or Nomoto et al. ( 120061 ), since the ta- 
bles cover only this range of elements. Element abundances for 
the Solar System are available from either spectroscopic abun- 
dance determinations from the solar photosphere or from labora- 
tory analysis of primitive meteorites. For the solar photosphere 
the table gives data from the compilation of Asplund, Grevesse 
& Sauval (2005 ). For meteoritic abundances the table shows data 
from the compilation of Palme & Jones (12003I) . 

The tabular values for the photosphere consider the recent 
significant downward revision of the abundances of O, C, and N 
by Allende Prieto, Lambert & Asplund (2001, 2002 ) compared 
to the previous compilations of Grevesse & Sauval (1998) and 
Anders & Grevesse ( 1989). The table also gives the abundances 
for C, N, and O derived by Holweger (2001 ), who also found a 
reduction of the solar abundances for these elements to be neces- 
sary, but not as much reduced as in the papers by Allende Prieto 
et al. The abundances of Allende Prieto et al. pose serious prob- 
lems to solar helioseismology (Delahaye & Pinsonneault 2006 ; 
Basu et al. 120071 ) while the abundances of Grevesse & Sauval 
( 1998 ) give good fits to observations. The incompatibility of the 
new C, N, O abundances with helioseismological results should 
be taken serious and the abundance reductions following from 



Table 2. Solar system element abundances a and the standard er- 
ror of the abundance determination cr. Abundances for the solar 
photosphere are from the compilation of Asplund et al. (2005 ), 
except for He where the recommended value for the early sun 
from Grevesse & Sauval d 19981 ) is given. Abundances for me- 
teorites are from the compilation of Palme & Jones (2003 ). For 
C, N and O also abundances from Holweger (2001 ) are given. 
For Ne and Ar see text. The last column indicates if solar pho- 
tospheric abundances from Asplund et al. (S), or from Holweger 
(H), or meteoritic (M) abundances are preferred for comparison 
with model calculation results; a small letter indicates that the 
element is not used for a comparison in the present work. 
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0.01 
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6.37 
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0.01 
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Si 
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19 


K 
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5.11 
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Ca 
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0.01 
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Ti 
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0.06 


4.95 


0.04 
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23 
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4.00 


0.02 


3.99 


0.02 
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24 


Cr 


5.64 


0.10 


5.67 


0.01 
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25 


Mn 


5.39 


0.03 


5.51 


0.01 
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26 


Fe 


7.45 


0.05 


7.49 


0.01 
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27 


Co 


4.92 


0.08 


4.90 


0.01 
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28 


Ni 


6.23 


0.04 


6.23 


0.02 


M 


29 


Cu 


4.21 


0.04 


4.28 


0.04 


M 


30 


Zn 


4.60 


0.03 


4.66 


0.04 


M 



using numerically calculated models for the solar convective 
flows (Asplund et al. [2000) to determine spectral line profiles 
seem to result in unrealistically small abundances. Also compar- 
ison with abundances in nearby G stars seem to indicate this (see 
section [2".3.5l ). The spatial resolution of their flow calculations of 
about 50 km (Asplund et al. 2000) compared to a pressure scale- 
height of the solar photosphere of about 125 km probably is in- 
sufficient and does not allow to account for velocity fluctuations 
on length scales small compared to the mean free path lengths 
of photons and therefore produces incorrect equivalent widths. 

For the solar photosphere no reliable abundances of the no- 
ble gases can be determined. For He a photospheric abundance is 
given in the table which is the value recommended by Grevesse 
& Sauval (1998 ) to be taken as the value of the He abundance 
of the early sun; the He abundance of the present sun is lower 
due to segregation effects and cannot be used for a comparison. 
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Fig. 7. Calculated element abundances relatively to solar abun- 
dances at the instant of Solar System formation (data according 
to Table[2]). Thin dotted lines show an as much as twice deviation 
from observed values 



The abundances for the other noble gases given in the table are 
determined from the Ne/Mg and Ar/Mg abundance ratios deter- 
mined from coronal lines as given by Feldman & Widing (2003 ). 
It is not sure that they really correspond to the initial solar abun- 
dances. 

For meteorites the abundance of H and of the noble gases do 
not reflect their abundance in the material out of which the par- 
ent bodies of the meteorites formed since these elements are not 
incorporated into the bodies of the early Solar System. Therefore 
no data for meteoritic abundances are given in the table for these 
elements. 

For meteorites the abundances of the volatile elements C, N, 
and O also are not representative for the abundances in the early 
Solar System since these elements are not (N) or only to a small 
fraction (C, O) condensed into solids and incorporated into the 
parent bodies of the meteorites. Correspondingly, the meteoritic 
abundances of C, N, and O given in the table are much lower 
than the photospheric abundances. For these elements the solar 
photospheric abundances have to be used for comparison pur- 
poses. A number of elements are highly volatile (cf., e.g., Palme 
& Jones 2003 ) and it is doubtful that these elements are com- 
pletely condensed in the parent bodies of the meteorites. Besides 
H, the noble gases, and C, N, and O, these are the elements CI, 
Br, I, In, Cs, Hg, Tl, Pb, Bi, from which CI is one of the elements 
in the table. For comparison purposes one should therefore use 
in principle the CI abundance from the photosphere, but since 
the abundance determination of CI for the solar photosphere is 
rather inaccurate, CI is presently not suited for comparison pur- 
poses. 

For the remaining elements, the meteoritic and the solar pho- 
tospheric abundances well agree, except for a number of heavier 
elements not contained in our table. For comparison with the re- 
sults of the chemical evolution calculation, we usually preferred 
the more accurate meteoritic abundance, whereas the solar pho- 
tospheric was used when both methods were of only moderate 
accuracy (as specified in last column of Tabled 

Additionally, the elements Li, Be, B, F are excluded from 
the comparison, since their production mechanisms are not im- 
plemented in the model program. 

In Fig.|7]we present the predicted element abundances of the 
ISM relative to Solar System abundances at the instant of Solar 
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Fig. 8. Calculated element abundances relative to abundances of 
F & G stars from the solar vicinity with with ages less than 1 Gyr 
(plusses), and of B stars from the range J"q ± 2 kpc of galactocen- 
tric distances (crosses) (data according to Table [3]). Thin dotted 
lines show an as much as twice deviation from observed values 

System formation at r — tq. Thin horizontal lines indicate a 
deviation by a factor of two upward or downward from Solar 
System abundances. As can be seen, the model fits the observed 
abundances with good accuracy. Most calculated element abun- 
dances reproduce the Solar System element abundances within 
a factor of about two, many elements even much better. The 
somewhat worse results for CI, K, and Sc have also be found 
by Kobayashi et al. (2006); the rather bad results for these ele- 
ments are of no importance for our model, since they are not one 
of the main dust forming elements. 

2.3.5. Abundances of F & G stars 

Abundances of the ISM cannot be measured directly because 
in the ISM the refractory elements are condensed into dust (cf. 
Savage & Sembach 1 19961 1. One possibility to determine indi- 
rectly total element abundances in the ISM is to determine at- 
mospheric abundances from some kind of 'young' stars which 
have not changed their surface abundances since they formed 
from interstellar matter. Best suited for this purpose are probably 
F and G main sequence stars from the galactic neighbourhood, 
that show the kinematics of thin disk stars and high metallicities, 
or, if stellar ages have been determined, are of an age of no more 
than a few Gyrs. For such stars one can assume that they formed 
from interstellar material with essentially the same properties as 
the present day ISM of the galactic neighbourhood. Bensby et al. 
(2005 ) and Bensby & Feltzing (2006) determined recently abun- 
dances for a number of elements for thin and thick disk stars 
from the solar neighbourhood. From the elements considered in 
that papers the following are relevant for our purposes: C, O, Na, 
Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, and Zn. 

First we consider from this sample the stars with a [Fe/H] 
ratio within +0.05 of the solar value. There are 6 stars which sat- 
isfy this condition and Table [3] shows the average abundances a 
of the above elements and the average scattering cr, of the abun- 
dances around the mean value. For comparison the table also 
shows the random errors of the abundance determinations from 
stellar spectra as given by the authors. These abundances are sur- 
prisingly close to the Solar System abundances, though the gen- 
eral metallicity is somewhat higher. If the range of metallicities 
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Table 3. Average abundances a of F and G stars with solar metallicity ( |A[Fe/H]| < 0.05 ) and of young stars (age < 1 Gyr) from 
the solar neighbourhood. <x abd is the accuracy of the abundance determinations from stellar spectra, o% is the scattering of the stellar 
abundances around the mean. Z is the metallicity calculated from these abundances 
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is increased to ±0.1 of the solar [Fe/H] ratio, the number of stars 
increases to 13, but the average values for the abundances are 
practically unchanged, i.e., the average abundances given for so- 
lar like stars in the table do not depend substantially on the pre- 
cise choice of the limit A[Fe/H]. Hence abundances of F & G 
stars with Solar System metallicity agree rather well with Solar 
System abundances as given in Table [2] that are already com- 
pared in Fig. |7]with results of our model calculation. Our model 
therefore reproduces reasonably well the observed abundances 
of solar metallicity stars at the solar cycle. 



Second we choose from the sample of Bensby et al. d20051 > 
and Bensby & Feltzing (2006 ) the thin disk stars with ages less 
than 1 Gyrs. Despite the large uncertainties of such age deter- 
minations it seems likely that these stars belong to the most re- 
cently born stars of the sample of thin disk stars. Their abun- 
dances therefore should sample the abundance of the ISM in the 
solar vicinity during the last, e.g., 1 ... 2 Gyrs or so. The av- 
erage abundances of the elements determined by Bensby et al. 
(2005 ) and Bensby & Feltzing (2006) for these stars are given 
in Table [3] Extending the sample to stars with an age less than 2 
Gyrs does not result in significant changes of the average abun- 
dances, i.e., the results do not depend on the precise choice of 
the age limit. The typical metallicity Z of the 'present' ISM de- 
termined in this way (the contribution of N and Ne to Z is es- 
timated) is slightly higher than the Solar System metallicity, as 
one may expect from ongoing element synthesis. Our model re- 
sults for the present day ISM abundances are compared in Fig. [8] 
with the observed abundances of F & G stars formed within the 
last Gyr given in Table[3] Our model results for the present ISM 
are also in good agreement with observations. 



2.3.6. Abundances of B dwarfs 

Stars of early spectral type B have short lifetimes. Therefore they 
sample abundances from the present day thin disk. Abundances 
have been determined in particular for B stars in stellar clus- 
ters and we show in Table [3] average abundances taken from the 
literature for B dwarfs in stellar cluster with galactocentric dis- 
tances from a range of ±2 kpc around the solar cycle. Despite 
the rather heterogeneous observational material the scattering of 
observed abundances around the mean is moderate, i.e., element 
abundances in the ring 8 + 2 kpc around the galactic centre seem 
to be quite homogeneous. The average abundances and, thus, 
the metallicity Z, are typically slightly less than the present-day 
abundances found from F and G stars (see Fig. [8]), as it is also 
found by Sofia & Meyer ( 12001 \ . Figure [8] compares our calcu- 
lated abundances for the present day ISM with the abundances 
of B dwarf; the agreement again is reasonable, but compared to 
the case of F & G stars it is worse since abundances of B dwarfs 
are lower than of F & G stars. 

2.3.7. Abundances of H II regions 

Abundance determinations from HII regions around massive 
stars give also some information on the element abundances of 
the present day interstellar medium, but only a small number of 
elements can be analysed. A serious problem is that the refrac- 
tory elements are in part or even completely condensed into dust 
particles, but the degree of depletion usually cannot be deter- 
mined with a sufficient degree of reliability, i.e., the total element 
abundance remains uncertain. Therefore we consider here only 
one case, the Orion nebula, since the abundances in this neb- 
ula should reflect element abundances in the present day ISM 
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Table 4. Dust species considered in the model calculation 
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Fig. 9. Characteristic cycle time of matter between the interstel- 
lar medium and stars at the solar cycle 



in the solar vicinity. Esteban et al. (2004) determined for the 
Orion nebula accurate element abundances for some elements. 
Table [3] gives their results for some elements which include al- 
ready their estimated corrections with respect to dust depletion. 
These abundances are almost the same as for young F & G stars. 
For the other refractory elements the degree of depletion is high 
(Shuping & Snow 1997 ) and they cannot be used for compari- 
son. Our model is also in accord with these abundances. 

The observational tests show that the rather simple type of 
model for the galactic chemical evolution considered here is al- 
ready suited to describe the global properties of the disk evolu- 
tion and in particular to describe the evolution of element abun- 
dances with sufficient accuracy that it may be taken as a basis 
for model calculations of the evolution of the dust content. 

2.4. Cycle time for the interstellar matter 

In the context of dust evolution an important quantity is the 
timescale for conversion of interstellar matter into stars. This is 
given by 

2ism 



^"cycle — 



B 



(8) 



This quantity is shown in Fig. [9] At the time of formation of the 
Solar System the cycle time of matter between the interstellar 
medium and stars was about 2.5 Gyrs. 

3. Production rates for Stardust 

In the following we study the evolution of the interstellar dust 
in the Milky Ways disk within the frame of the simple one- 
zone approximation for galactic disk evolution described above. 
The model follows the principles of the model of Dwek ( 1998 ), 
which was the first one that coupled dust evolution consistently 
with a full model for the chemical evolution of the Milky Way. 
The present model calculation concentrates on the most abun- 
dant dust components formed in stars that are found also as 
presolar dust grains in meteorites: 

- silicate dust 

- carbon dust 

- silicon carbide. 

It also considers 



- solid iron 

though that is not yet identified as a presolar dust species. 
For theoretical reasons it should, however, be an abundant dust 
species produced in stellar outflows. 

3.1. Evolution model for the interstellar dust 

In our evolution model we consider a number of different dust 
species, denoted by an index j. 

First we discriminate between dust coming from different 
types of parent stars. Even if the chemical composition of a cer- 
tain dust species formed in outflows or ejecta of different stellar 
types is the same, the individual grains of this dust species are 
carriers of the isotopic anomalies corresponding to the particu- 
lar nuclear processes operating in their parent stars. If they are 
investigated in the laboratory as presolar dust grains, one can, at 
least in principle, identify for every grain its formation site. This 
makes it desirable to count dust species from different types of 
stars with the same overall chemical composition but with differ- 
ent types of isotopic anomalies as different species j. Presently 
we choose a not too fine division into stellar types and consider 
three different kinds of stellar sources: Supernovae of type II, 
supernovae of type la, and AGB stars. Supernovae of type II and 
AGB stars can form all of the four chemically different types of 
dust considered in our model, while supernovae of type la prob- 
ably can form only iron dust (if ever). Hence we consider in our 
model nine different kinds of Stardust coming from three differ- 
ent types of parent stars. 

Second we consider dust formed in the interstellar medium 
itself. From the element abundances in the interstellar medium 
one expects that silicate dust can be formed. Observations of the 
interstellar dust indicate, that also carbon dust can be formed 
in certain regions of the interstellar medium. It seems unlikely 
however, that SiC dust can be formed, since this requires a car- 
bon rich environment, which is not encountered in interstellar 
space. Iron dust may also be formed in the ISM, though this ele- 
ment is probably mainly consumed by silicate formation. Hence 
we consider in our model three kinds of dust formed in the inter- 
stellar medium: silicate, iron, and carbon dust. 

Totally, our model considers the twelve different kinds of 
dust from stellar sources and the interstellar medium given in 
Table H 

We describe the abundance of each dust component j in the 
interstellar medium by its surface mass density The evolu- 
tion of the surface density is determined by the equation 



dE 



7,d 



dr 



-ISM 



B 



+ ^ RtfA - L J,d + Gj, d . 



(9) 



The first term on the rhs. describes the loss of dust from the 
interstellar medium by star formation. It is assumed that only the 
dust content of the matter that is converted into stars is destroyed 
and that no additional dust is destroyed during the course of this 
process. Also no return of freshly formed dust from protostellar 
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disks by winds or jets is assumed to occur, though this has been 
speculated to be possibly important (Tielens 2003). The second 
term on the rhs. describes the amount of dust of kind j injected 
by stars of type I into the interstellar medium. The third term 
on the rhs. describes the losses, the destruction of dust of kind j 
in the interstellar medium, mainly by supernova shocks, and the 
last term describes the gain, the formation of dust of kind j in the 
interstellar medium by growth processes in molecular clouds. 

From equations (0 one calculates the surface mass density 
of the different dust species j. Additionally one has the set of 
equations for the total surface densities S, of each element i in 
the ISM. The surface density E; >g of each element in the gas 
phase of the ISM then follows as the difference between its total 
surface density X, and the sum of the contributions of all dust 
species containing that element 

stg = s*-E v yr Lz ^- (10) 

Here vy is the number of atoms of element i in one formula unifl 
of dust species j, Ai and Ajj are the atomic weights of element 
i and of one formula unit of dust species j, respectively, and the 
summation is over all dust species containing element i. 

In the following we describe the details of the processes rel- 
evant for the evolution of interstellar dust abundances and how 
they are implemented in the model. 

3.2. Dust return by AGB stars 

The main factories of dust in space are low and intermedi- 
ate mass stars in the AGB stage of their evolution. These are 
stars with initial masses between about 0.8 Mq and about 8 Mq, 
which end their life as White Dwarfs. The lower mass limit cor- 
responds to that initial mass, for which the lifetime of a star cor- 
responds to the age of the Milky Way. The upper mass limit cor- 
responds to stars which finally explode as supernovae and do not 
evolve through an AGB phasfl 

The initial element mixture of all stars is oxygen rich in 
the sense that the abundance eo of oxygen exceeds the abun- 
dance ec of carbon. This doesn't change during their whole 
evolution up to the thermally pulsing AGB (TP-AGB), despite 
some abundance changes during first and second dredge-up on 
the Red Giant Branch and the early AGB, respectively. If the 
third dredge-up process starts operation, the ashes of He burn- 
ing are mixed to the convective envelope of the star after each 
thermal pulse, increasing stepwise the carbon abundance of the 
convective envelope, but only marginally changes its oxygen 
abundance. The resulting evolution of the carbon to oxygen ratio 
fc/^o on the TP-AGB depends on the initial mass of the stars: 

(1) Low mass stars with initial masses less than about M = 
1.5 M© loose their envelope by a stellar wind before the car- 
bon abundance exceeds the oxygen abundance. Only for low 
metallicities their envelope becomes carbon rich prior to com- 
plete envelope loss. For very small initial masses probably no 
third dredge-up occurs and the stars never become carbon stars. 

(2) Stars from the range of initial masses between about M = 
1 .5 M© and about M = 4 M© increase after a number of thermal 
pulses their carbon abundance over the oxygen abundance and 

2 The formula unit is the fictitious molecular group in the solid cor- 
responding to the chemical formula of the condensed phase 

3 Presently we neglect the possibility that stars from the region of 
initial masses 8 < M < 12Mq may become super- AGB stars (e.g. 
Gil-Pons & Garcfa-Berro l2002) > 



become carbon rich. They further evolve as carbon stars until 
their envelope is lost by the stellar wind. 

(3) Intermediate mass stars with initial masses between about 
M = 4 M© and about M — 8 M© convert the dredged-up carbon 
rapidly into 14 N via the CN-cycle, since the lower part of their 
convective envelope overlaps with the upper edge of the H burn- 
ing shell. The oxygen, however, is not affected by this process, 
and for this, the carbon abundance in the envelope of these stars 
is much less than the oxygen abundance. They do not become 
carbon stars until most of their envelope is lost and the convec- 
tion zone of the remaining envelope no longer overlaps with the 
H burning shell. For a short period the stars then also become 
carbon stars until finally the last portion of their envelope is re- 
moved by the stellar wind. 

Depending on the C/O abundance ratio, the stars produce 
different dust mixtures in their outflows. Here we are only inter- 
ested in the dominating dust species, which are formed from the 
most abundant elements and concentrate on the following four 
types of dust: silicates, carbon, silicon carbide, and iron. The 
dust production rate over the whole period of AGB evolution for 
these dust species was calculated by Ferrarotti & Gail (120061 ) . 

The silicate dust is only produced during the oxygen rich 
phase of the stellar evolution where the stars spectroscopically 
appear as M stars. Some minor fractions are also produced dur- 
ing the S star phase where the C/O abundance ratio is close to 
unity. The silicates are a mixture of olivine- and pyroxene-type 
amorphous dust and for part of the stars also up to about 15% 
of nearly iron free crystalline forsterite and enstatite is observed 
to be formed (cf. the review of Molster & Waters 2003). The 
present work does not discriminate between the two types of 
amorphous silicate dust since for silicate dust in the interstel- 
lar medium it is presently not possible to discriminate by ob- 
servations of the dust absorption spectrum unambiguously be- 
tween the two different components (cf. the contradictory results 
in Chiar & Tielens 120061 and Min et al. 120071 ). Also crystalline 
Mg-silicates are not considered since they are not found in the 
interstellar medium (Kemper, Vriend & Tielens 2004 ), possibly 
because they are rapidly amorphized in the ISM after their ejec- 
tion by interaction with energetic electrons and ions (cf. Demyk 
et al. 120041 Jager et al. 120031 . 

Carbon and silicon carbide dust are produced by AGB stars 
during their carbon rich phase of evolution on the AGB where 
they spectroscopically appear as C stars. 

Iron dust is included in the model calculation, though it has 
not yet unambiguously been identified as major dust component 
in stellar outflows; only some hints for its existence have been 
found up to now (e.g. Kemper et al. [2002). This is because no 
readily identifiable spectroscopic features exist for solid iron. 
Nevertheless, for reasons of element abundances it should be 
an abundantly formed dust species in S stars and C stars and 
to some extent also in M stars. 

MgS is also observed to be an abundant dust component in 
many C stars (cf. Molster & Waters 2003 ), but is not included 
in the modelling, since it is not yet clear by which mechanism 
it can be formed in a stellar outflow. This, of cause, presently 
prevents its modelling. 

Figure [10] shows the calculated dust production rates for 
the four types of dust considered. In the model calculation of 
Ferrarotti & Gail (2006), olivine- and pyroxene-type dust are 
treated as separate species but their production rates are added 
for the reasons mentioned above. There is a general tendency for 
the stars either to be a factory mainly for silicate dust or to be a 
factory mainly for carbon dust (cf. Fig. 12 of Ferrarotti & Gail 
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(a) 






Fig. 10. Dependence of the dust masses returned by single AGB stars for the four main kinds of dust species (silicates, carbon, 
silicon carbide, and iron) on metallicity Z and initial stellar mass M*. All masses are in units of Mq. Data from Ferrarotti & Gail 
(2006 ) with some additional models 



2006), because most of the dust formed over the total lifetime of 
a star on the AGB is formed during the very last pulse cycles on 
the TP-AGB, where mass-loss rates are highest. If the stars are 
carbon stars during this phase they mainly produce carbon dust 
(and SiC), otherwise they mainly produce silicate dust. 

The carbon dust production, shown in Fig.fTOh. is dominated 
by stars with initial masses between about 1.5 and 4Mq and 
does not vary much with initial stellar metallicity, since the car- 
bon required for carbon dust production is synthesized from He 
by the star itself. Stars with initial masses M > 4 M0 also form 
some carbon dust, but only small amounts during their very last 
stage of evolution when hot bottom burning is no more active. 
Stars with initial masses M < 1.5 M0 do not form much carbon 
dust because the total mass returned by them on the AGB is quite 
small or because they do not suffer third dredge-up events (for 
very low initial masses) or a too small number of them. 

The production of the other dust species by AGB stars 
strongly depends on their initial metallicity because the required 
heavy elements — with the possible exception of Mg — are not 
fabricated by AGB stars but have to be formed in many preced- 
ing stellar generations until their abundances grew to a level that 
dust formation becomes possible. 



Figure [TOb shows the silicate dust production by AGB stars. 
The silicate production is efficient for stars from essentially that 
range of initial masses where they do not become efficient car- 
bon dust producers, i.e., the main contribution comes from stars 
with initial masses M < 1.5 M0 or M > 4M . But also in the 
mass range in between, where the stars are efficient carbon dust 
factories, they produce some silicate dust before they become 
carbon stars. The silicate dust production starts to become effi- 
cient only at rather high metallicities because only then sufficient 
amounts of Si, Mg, and Fe for silicate formation are available in 
the stellar outflows. 



Figure [TOfc shows the production of silicon carbide dust by 
AGB stars. This is produced by carbon stars and therefore its 
production is limited to the same range of initial masses as for 
carbon production. The lack of available Si, however, also pre- 
vents the formation of much SiC in low metallicity stars. 

Figure [TOtl shows the production of iron dust by AGB stars. 
Iron dust formation seems to be efficient in outflows from AGB 
stars only at rather high metallicities which are not encountered 
in the Milky Way at the solar cycle, but only close to its centre. 
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The dust-mass injection rate of dust species j into the inter- 
stellar medium is given by 

Rj,i,i(t,r) = dM <&(M) —tj— 

Jm, Mm 

M m {M,Zi$ u {t h ,r)). (11) 

The index I here refers to AGB stars. The quantity Q> is the initial 
mass function, described in Sect. 12.1.41 M av is the average mass 
of the stars, given by Eq. ©. B{t, r) and Zism(?, r) are the stellar 
birthrate and the metallicity of the interstellar medium at instant 
t and galactocentric radius r, respectively. Both quantities are 
taken from the model calculation for the evolution of the galactic 
disk (see Figs.[2]and[3]l. The instant 

t b = t-T(M,Z 1SM (t h ,r)) (12) 

is the time of birth of a star ending its life at instant t . The quan- 
tity t(M, Z) is the lifetime of a star with initial mass M and 
metallicity Z. The metallicity of the stars equals the metallic- 
ity Zism of the interstellar medium at their birthtime fb, that is 
taken from our model for the Milky Ways evolution (cf. Fig. [3]). 
Note that % is given by a non-linear equation which has to be 
solved numerically for each M and t. Tables for the mass return 
in the dust species j by AGB stars of different initial masses and 
metallicities are given in Ferrarotti & Gail (2006); some addi- 
tional models have been calculated for the present work. The 
integration is performed over the initial masses of the stars from 
the lower limit M\, here taken to be 1 Mq, and the upper limit 
Mwd up to which stars evolve into White Dwarfs, here taken to 
be8M Q . 

3.3. Dust production by massive stars 

In principle one has three different processes contributing to the 
dust return by massive stars that finally explode as SNe: 

1. Dust formed in the massive cool stellar winds of Red 
Supergiants, i.e., massive stars on the Red Giant Branch. 
This is relevant only for stars from the range of initial masses 
8 ^ M % 40 Mq, since only stars from this mass range enter 
the Red Giant stage. 

2. Dust formed in massive shells of 1 ... 10 Mq ejected by re- 
peated giant eruptions during an LBV-phase, like for in- 
stance that observed in 77 Car (cf. Smith & Owocki 2006). 
This is relevant only for massive stars from the region of ini- 
tial masses M <: 40 Mq which evolve through a LBV phase. 
Most of the mass ejected by these very massive stars prior to 
their SN explosion seems to be ejected in a few such events 
(Smith 120061 1 which are accompanied by copious dust for- 
mation. 

3. Dust formed in the ejected matter after the final supernova 
explosion. 

The dust grains formed in outflows from Red Supergiants and 
giant eruptions of very massive stars carry isotopic anomalies 
resulting from hydrogen burning via the CNO-cycle, while the 
dust grains formed after a supernova explosion show the very 
different isotopic signatures from heavy element synthesis. Both 
types would be clearly distinguishable, if investigated as preso- 
lar grains in the laboratory. Both types should be included as 
separate types of dust in a model calculation. 



3.3.1. Dust formed in pre-supernovae 

The dust formed in stellar winds or ejecta prior to the supernova 
explosion is later subjected to the shock wave from the SN ex- 
plosion. This shock wave destroys the dust in the swept-up ma- 
terial if the expansion velocity exceeds 150 kms~' (e.g. Jones 
et al. |19961 l. For a simple estimation of the importance of this 
process we consider the case that the blast wave expands into a 
medium with constant density. At the end of the adiabatic ex- 
pansion phase the radius and the velocity of the shocked region 
are about (e.g. Shull & Draine [T98"7j l 

R sh = 16.2 E 2 5 [ 7 n vi pc (13) 
V sh = 331 El[ H n n km s- 1 . (14) 

The transition between the Sedov-Taylor expansion and the sub- 
sequent pressure dominated snowplow phase occurs at 

t S T-PDE = 1-91 x 10 4 E 3 5 [ H n Q 4n yr. (15) 

Here £51 is the explosion energy in units of 10 51 erg and «o the 
density of the ambient medium in units of 1 cm 4 . Since in the 
snowplow phase the shock velocity drops rapidly, the dust de- 
struction occurs mainly up to the end of the adiabatic expansion 
phase given by Eq. HS\ . 

Red Supergiants: First we consider the case of Red 
Supergiants and let £51 = no = 1. Typical expansion veloci- 
ties of stellar winds of supergiants are v exp = 20 kms -1 . The 
wind material requires a time of about R s h/v eX p - 790000 yrs 
to expand to the distance R^. The shock strength then is suf- 
ficient to destroy all dust material ejected during a period of 
7.9 x 10 5 - t ST -pde ~ 7.7 x 10 yrs before the SN explosion. 
This is of the order of the evolution time on the Red Giant branch 
(e.g. Schaller et al. 1992). The main period of dust formation of 
such stars, however, is much shorter. Mass-loss rates of super- 
giants during the phase where they are enshrouded by massive 
dust shells are of the order of 10~ 4 . . . 10~ 3 Mq yr -1 (e.g. van 
Loon et al. 1 19991 ) and this phase can last at most about 10 5 yrs 
otherwise the stellar envelope over the He core would be lost 
completely by the stellar wind prior to explosion, which is not 
observed for this mass-range. 

Hence, all dust formed by Red Supergiants is expected to be 
destroyed by the shock wave of the subsequent supernova explo- 
sion. Even if some dust survives in some cases, Red Supergiants 
cannot be important sources for interstellar dust. 

Luminous Blue Variables: The expansion velocity of the 
matter from giant eruptions is somewhat higher than for winds 
of Red Supergiants (cf. Lamers et al. 2001) and may be as high 
as 100 kms" 1 . Correspondingly, the supernova shock destroys 
all the dust that was ejected by a giant eruption if the supernova 
explosion follows within about 2 x 10 5 yrs after the end of the 
LBV phase. The LBV phase, however, seems to occur during 
the first transition from the blue to the red in the Hertzsprung- 
Russel diagram (Lamers et al. 2001) and is followed by a WR- 
phase which lasts about 3 x 10 5 . . . 10 6 yrs (Meynet & Maeder 
|20051 l. If the star finally explodes, the velocity of the shock from 
the SN explosion is already too slow for destroying the dust at 
the instant when it catches up with the ejected LBV shell. Dust 
formed in giant eruptions could therefore be an important source 
for interstellar dust. Unfortunately, however, there is presently 
insufficient information on the dust production by these objects 
to include them in the model calculation and dust production by 
LBVs therefore cannot be considered in our present model cal- 
culation. 
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Clearly, the real situation is more complex since a supernova 
explodes into the matter ejected by the stellar winds of the pre- 
ceding evolutionary stages (cf. Dwarkadas 2006 for a brief dis- 
cussion), or into the hot bubbles of other supernovae, and the 
dust ejected by one massive star may be subjected to the SN 
blast waves of other massive stars from the same stellar cluster. 
A more detailed investigation of the whole problem is required 
to determine the survival probability of dust formed by a star 
prior to its SN explosion. 

3.4. Dust return by supernovae 

Unfortunately it is presently not definitely known which super- 
novae do form dust and in which quantities. Undoubtedly there 
is some dust formed by supernovae since presolar dust grains are 
known which bear the signatures of element synthesis in super- 
novae. The abundance of X-grains in the population of preso- 
lar SiC grains, however, is small compared to mainstream SiC 
grains (cf. Hoppe et al. 2000; Nittler & Alexander [2Q031 > which 
are thought to come from AGB stars. Dust formation by super- 
novae, therefore, seems to be an inefficient process. For super- 
novae of type la, observations even seem to indicate that they 
do not form dust at all (Borkowski et al. 120061 ). From the theo- 
retical side also little is known about dust condensation in SNe; 
only a few model calculations for dust condensation in super- 
nova ejecta are available (Kozasa et al. 119891 Todini & Ferrara 
I200T1 Nozawa et al. 120031 Schneider et al. 120041 1 and they are of 
a very qualitative nature. 

Presently there are no reliable models for dust formation in 
supernovae available on which one can base a modelling of the 
contribution of supernovae to the interstellar and presolar dust 
population. Therefore we apply in the present model calcula- 
tion the same simplified approach as in Dwek (1998 ) to account 
for the contribution of supernovae to the dust production in the 
Milky Way. It is assumed that supernovae of type II produce all 
types of dust considered here and supernovae of type la produce 
only iron dust. The dust return rate is assumed to be given by 
the total mass return rate of the key element required to form 
a particular kind of dusQ times some efficiency factor rj. This 
efficiency factor is simply guessed or is estimated from observa- 
tional quantities. 

We therefore use the following production rates for the dust 
species 



sil 



fisil,/,d(f, >") = ?7sil,SNII^Si,SNIl(f,'") — 

Asi 

A C ar 

^car,/,d(f, r) = ?7car,SNII^C,SNn(*,'') — — 

Ac 



^sic,/,d(f> r ) - 77sic,SNII^Si,SNIl(f. r ) 



A sic 

A- 

Riro.lA 1 ' r ) ~ / 7ii-o,SNII^Fe,SNIl( f >'") ~. ■ 

AFe 



(16) 
(17) 
(18) 
(19) 



The index I here refers to supernovae of type II. ^si.SNii is the 
rate of mass return to the interstellar medium of element Si by 



4 As key element we usually choose that one of all the elements 
forming the considered dust species, for which the quantity e/i takes 
the smallest value. Here e is the abundance of an element by number 
relative to H, and ( is the number of atoms of the element in the chemical 
formula of the compound. The key element determines the maximum 
amount of dust material that can be formed for the considered species 



Table 5. Characteristic quantities and numerical coefficients 
used for calculation of grain destruction and grain growth (SiC 
does not form in the ISM) 





silicates 


carbon 


iron 


SiC 


T;,snr (Gyr) 


0.4 


0.6 


0.6 


0.6 


^SNII 


0.00035 


0.15 


0.001 


0.0003 


T/SNIa 


0.0 


0.0 


0.005 


0.0 


key element 


Si or Mg 


C 


Fe 


Si 


Atomic weight A 


121.41 


12.01 


55.85 




bulk density p c (gem -3 ) 


3.13 


2.25 


7.86 





all supernovae of type II, defined by 

B(t b ,r) 



flsi,SNii(f, r) 



I ' 

JM\ VD 



dM<S){M)- 



Msi, re t(M,Zis M (fb,r)). 



(20) 



The quantity Msi, re t(M, Z) is the Si mass returned by a supernova 
with initial mass M and metallicity Z. The integration is from the 
lower limit Mwd to the upper limit M u , here taken to be 40 Mq 
since the tables for mass-return by supernovae of Woosley & 
Weaver ( 1995 ) and Nomoto et al. (2006) only extend up to this 
upper mass. The mass return rates for the other elements are 
defined correspondingly. 

The quantities A s n, A car , A s ; c , A[ m are the atomic weights of 
the dust species and A$i, A car , Ap e the atomic weights of the key 
elements. 

For supernovae of type la the dust production rate is 



Rho,l,d(t, r) — 77i ro ,SNIa ^Fe,SNIa(f, r) — , 

AFe 



(21) 



where the index / now refers to supernovae of type la. The mass 
return rate for iron is 



#Fe,SNIa(*, r) = Af Fe , ret /?SNI a (*, r) , 



(22) 



where ^sNia(?, r) is the supernova rate of type la, which is taken 
from the galactic evolution model (see Fig.|2|. Mpe.ret is the iron 
mass returned by supernovae of type la and is taken from the 
tables of Iwamoto d!9991 >. 

The quantities /7 s fl,SNn, ^iro.SNia are the efficiencies for 
conversion of the key elements of the different dust species into 
dust particles and refer to the amount of dust injected into the in- 
terstellar medium in relation to the total mass of the key element 
returned to the interstellar medium. The dust first formed in the 
expanding SN ejecta is later overrun by the reverse shock and 
part of it is destroyed again (Dwek 120051 ). Further, part of the 
surviving dust may be later destroyed by blast waves from other 
supernova explosions in the same stellar cluster. The efficiencies 
77 as they are defined here consider all destruction effects and 
may therefore be significantly smaller than the efficiency of the 
initial dust condensation. 

So far only a few attempts have been made to estimate the 
condensation efficiency in SNe by analysing spectroscopic data, 
resulting in very different dust yields, from only 5 X 10~ 4 to 
4 x 10~ 3 for type II SN 1987A (Erco lano et al. [2007b to 0.12 
for SN 2003gd (Sugerman et al. 120061 The values of 77 for dif- 
ferent types of SNe are still unknown and have to be guessed 
somehow. The numerical values chosen in this paper are much 
smaller than the values assumed in Dwek (1998) and are given 
in Table [5] The values for the efficiencies ?7sU,snii> ?7sic,SNii, and 
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Fig. 11. Evolution of the dust injection rates at the solar cycle 
from different stellar sources 



^car.SNii of silicate, SiC, and carbon dust formation in SN II, re- 
spectively, are estimated from the abundances of presolar sili- 
cate grains from supernovae, of X-type SiC grains, and graphite 
grains from supernovae. This is discussed in Sect. 15.31 The ef- 
ficiencies 77i r0; sNii an d ^iio.SNia f° r i ron dust production in SNe 
of type II and type la, respectively, are arbitrarily set to a low 
non-zero value, but they may well be equal to zero. 

3.5. Dust injection rates for the Milky Way model 

Figure [TT] shows the variation with time of the dust injection 
rates from stellar sources into the interstellar medium at the solar 
cycle, as calculated for our model of the evolution of the Milky 
Way. The dust injection rate in this model is dominated by car- 
bon dust from AGB stars and SNe, and by silicate dust from 
AGB stars, except for the very first period before the first ap- 
pearance of AGB stars, where dust return from SNe dominates. 
The SN injection rates are very uncertain, however, since they 
depend on the efficiencies r\ which are only badly known and in 
this paper are determined from abundance ratios of presolar dust 
grains from AGB stars and supernovae. 



4. Dust evolution in the ISM 

4. 1 . Dust destruction in the interstellar medium 

Once condensed in a circumstellar envelope dust grains are 
driven away from the star by radiation pressure and are grad- 
ually mixed with the ambient ISM. The analysis of astronomical 
observations and presolar dust grains from meteorites indicate a 
high degree of subsequent dust processing in the ISM. The pro- 
cesses can be divided in two groups: 

1. the destruction processes (thermal sputtering, evaporation 
in high-velocity grain-grain collisions, chemical sputtering), 
and the growth by accretion that changes the total dust mass; 

2. the processes modifying the grain size distribution or the 
phase of the grains (shattering in grain-grain collisions, co- 
agulation, amorphization of crystalline Stardust etc.). 

Detailed discussions of dust processing in the ISM can be found 
in Tielens (2005), for example. Since in this paper we do not 
study the grain size distribution, we are interested in the first 



group of processes only, i.e., in the destruction processes and the 
growth processes in the ISM. 

At the stage of grain ejection from the star by radiation pres- 
sure, grain destruction is not important for quiescent mass loss. 
Typically, drift velocities for sources with high mass loss are 
about 1 km s _1 (e.g. Gail & Sedlmayr 1985), which is less than 
the threshold for sputtering of the dust; grain erosion is there- 
fore unimportant during dust ejection into the ISM (Woitke et al. 
1993 Jones et al. 2005). However, the lower limit for shattering 
in grain-grain collisions is about 1-2 km s _1 and grain fragmen- 
tation could be important if there are significant relative veloc- 
ities between the grains in stellar envelopes (Jones et al. 1996). 
In case of dust condensation in SN explosions freshly formed 
grains later are overrun by the reverse shock with velocities of 
the order of at least 1000 km s~ . Therefore a significant part of 
the SN condensates is destroyed by sputtering in reverse shocks 
(Dwek 2005). This is included in the efficiencies 77 introduced in 
Sect. El 

The dominant dust destruction process is dust destruction 
in the ISM by sputtering in high-velocity SN shocks (v > 
150 km s -1 ) (cf. Jones et al. 1 19961 1 that results in injection of 
atoms into the gas phase by interaction with impinging energetic 
ions. This process works almost exclusively in the warm phase 
of the interstellar medium, which links the dust destruction prob- 
lem in the ISM in principle closely with the multiphase structure 
of the ISM. Since we approximate in our present model the ISM 
by a simple one-phase model, we describe this process in the ap- 
proximation which describes the destruction process in terms of 
grain lifetimes against destruction by SN remnants t^snr, i-e., 
by the average time required to return the grain mass back to 
the gas-component of the ISM. The change of surface density of 
dust species of kind j per unit time by dust destruction is: 



-J4 



T/,SNR 



(23) 



Since we consider only equations for the total element abun- 
dances (gas+dust) in the ISM and no separate equations for the 
gas-phase abundance, the mass return to the gas phase needs no 
special treatment. This is automatically accounted for if the gas 
phase abundances are calculated from Eqs. (TTOb . 

The efficiency of the dust destruction process is determined 
by the SN II and SN la frequency, the density of the ISM, as well 
as by the intrinsic properties of grains such as size, composition, 
and bond energies. The grain lifetime at any radius r and instant 
t can be expressed by the current lifetime at the solar circle ac- 
cording to Dwek ( 1998) by the approximate scaling law: 



T/,SNR0,r) 



2iSM(fG,>"o) Rsn(t,r) 



(24) 



The total SN rate RsN(t, r) can be replaced by the stellar birthrate 
B(t, r) with the assumption of spatial and temporal constancy of 
the IMF. The dust species considered in our model are listed in 
Table H For these dust species we need to know the destruction 
timescales t,;snr(/, r). 

Grain lifetimes t^snr(^g> r©) have been subject of many the- 
oretical studies considering thermal sputtering, shattering and 
models of SN explosions; the most detailed ones are Jones et 
al. ([1994 1996 see also the references therein). For the current 
grain lifetimes at the solar cycle against destruction in all phases 
of the ISM we take the theoretical estimates of 0.6 and 0.4 Gyr 
for carbonaceous and silicate dust, respectively, from Jones et 
al. ([1996). Unfortunately they did not present the corresponding 
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timescales for iron and silicon carbide dust, but their results for 
the sputtered dust mass fraction for different shock velocities for 
iron and for silicon carbide dust are somewhat higher but similar 
to carbon dust. We therefore choose for both the same lifetime 
of 0.6 Gyr as for carbon dust. The model results depend only 
moderately on the value of t ^snr as long as t ; snr is at least of 
the order of 0.5 Gyr. This is because for this and higher values of 
T/ ; snr the destruction timescale for dust becomes roughly com- 
parable with the timescale of conversion of ISM matter into stars 
that is of the order of 2.5 Gyrs (cf. Fig. [9j». In the opposite case, 
for timescales t^snr significantly smaller than about 0.5 Gyrs, 
the results depend markedly on the precise value of t^snr- 

The values of t^snr(^g> r o) use d m the model calculation are 
also shown in Table[5] 

4.1 .1 . Evolution of Stardust 

Stardust is only destroyed in the ISM and does not gain mass by 
accretion of gas phase material. All material from such grains 
ejected into the gas phase rapidly mixes with the existing ISM 
gas-phase material and the specific isotopic anomalies carried 
by the Stardust material are lost by mixing together eroded ma- 
terial of grains from many different kinds of stellar sources. If 
such material is later accreted by dust grains in the ISM it shows 
no isotopic anomalies. Even if such material grows as mantle 
material on Stardust cores the differences in isotopic composi- 
tion between core and mantle survive since dust grains in the 
ISM are not expected to ever become hot enough (> 1000 K) 
for sufficient long periods that solid state diffusion smoothes out 
isotopic abundance differences between a Stardust core and an 
ISM-grown mantle. Hence any accreted mantle material can be 
clearly discriminated (if it could be analyzed in the laboratory) 
from cores originating from stellar sources by showing isotopic 
abundances ratios close to Solar System isotopic abundance ra- 
tios, even if the general chemical composition and mineralogical 
structure of an ISM-grown mantle material should resemble that 
of a core with stellar origin. Therefore we treat the dust species 
from stellar sources in our model as separate dust components 
and omit for these species in Eq. © the growth term; only the 
destruction and source terms are retained. 

4.2. Dust growth in the interstellar medium 

Dust grains cycle between the cloud and intercloud phase of the 
ISM on a timescale of 3 x 10 7 yr (e.g. Draine [T9901 Tielens 
1998 ), undergoing destruction in the warm intercloud medium. 
All theoretical calculations of grain lifetimes against destruc- 
tion by SN shocks agree that they are much shorter than the 
~2 Gyr timescale of dust injection by stars (e.g. Jones et al. 1996; 
Tielens et al. 2005 ). This requires an efficient mechanism of re- 
plenishment of the dust content of the ISM. Another evidence 
of dust growth in the ISM is, that gas abundances in the ISM of 
major dust forming elements show strong depletion in compari- 
son to solar abundances (e.g. Savage & Sembach |19961 Jenkins 
2004 ), which correlates with the ISM density. The only possible 
site of grain growth in the ISM are the dense molecular clouds 
of the cold phase of the ISM (Draine [1990] ). 

It is known that the density of the ISM is not sufficient to 
allow for formation of new dust grains, only low temperature 
accretion of refractory material on pre-existing stellar grains is 
possible. The thin mantles accreted in the ISM are likely more 
volatile than stellar dust and can be lost more readily during dust 
cycling between ISM phases. Besides, in dense molecular clouds 



the accretion will be faster, and grains will probably be formed 
far from equilibrium, therefore one would expect the grain man- 
tles to be amorphous and heterogeneous (Jones 2005 ). Thus, dust 
accreted in molecular clouds (the MC-grown dust) has proper- 
ties different from that of Stardust and is treated in our model as 
a separate component denoted by an index ISM. 

Dust growth in molecular clouds by accretion on existing 
grains needs to be considered in our model for silicate and car- 
bon dust, since it is known that observed depletions of the ele- 
ments in the ISM are too high as that only destruction can be ac- 
tive (cf. Jenkins 2004]for a recent discussion). The grains which 
serve as growth centres for accretion of gas-phase material need 
not necessarily be the Stardust particles, though these are needed 
to serve as initial growth centres for a start-up of the whole 
process. Also fragments formed from shattering of MC-grown 
grains by SN shock waves in the warm component of the ISM 
may serve as growth centres for accretion of refractory elements 
in the gas phase if mixed into molecular cloud cores. 

An unclear case for growth in the ISM is iron dust which 
might be a component of the ISM dust mixture if not all Fe is 
used up by the formation of magnesium-iron-silicates, in which 
case it also could grow in the ISM. However, metallic iron is 
probably unstable against oxidation in the ISM (Jones 1990), 
while, on the other hand, iron oxides seem not to form a signifi- 
cant species in the ISM dust mixture (Chiar & Tielens 2006}. We 
consider in our model calculation iron as a possible MC-grown 
dust component since there are observational indications that not 
all condensed iron always resides in silicates (cf. Cartledge et al. 
[20061 their Fig. 10). 

How the growth process works in detail is not definitely 
known. For interstellar carbon dust it may proceed in the way 
described in Jenniskens et al. (1993) as a multistep process, 
initiated by deposition of ice mantles, and proceeding via car- 
bonization and polymerization driven by UV irradiation. The 
problem of growth of silicate dust in the ISM has long time 
remained unsolved, because the formation of tetrahedral S1O4- 
structures probably requires higher temperatures, than the 10- 
30 K observed in molecular clouds. At these low temperatures 
ice mantles are likely formed on the grain surface, preventing 
further growth of silicates. The solution of the problem is possi- 
bly provided by intermittent dissipation of turbulence in molecu- 
lar clouds (Falgarone et al 2006 ). The large local release of non- 
thermal energy in the gas by short bursts of turbulent dissipa- 
tion has been shown to be able to trigger a specific warm chem- 
istry, which can be traced by the high abundances of CH + , H2O, 
and HCO + observed in diffuse gas. It is shown that signatures of 
warm chemistry survive in the gas more then 10 3 yr during chem- 
ical and thermal relaxation phase, see Fig. 10 in Falgarone et al. 
(2006). Such local change of the gas temperature could provide 
the mechanism for further silicate growth, if the grain tempera- 
ture increase is sufficient for ice mantles to evaporate. The latter 
is defined by equating energy from collisions with warm gas and 
the emitted infrared energy, and thus depends on infrared absorp- 
tion coefficients of mantle and core grain material, e.g. it differs 
noticeably for water and organic ices. 

Our preliminary estimates show that the energy released lo- 
cally by turbulent dissipation in molecular clouds is sufficient to 
evaporate organic ice mantles from the surface of silicate grains, 
although detail calculations of temperatures and residence time 
in the relaxation phase for grains with different compositions 
have to be done to make quantitative estimates. This is a sep- 
arate problem important for understanding the physics of dust 
growth in molecular clouds will be studied in further papers. 
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4.2.1 . Growth of dust grains in molecular clouds 

In calculating the growth rates for the dust species we follow a 
different procedure as in Dwek ( 1998). Some modifications are 
necessary because (1) we wish to consider specific dust com- 
ponents and not merely the surface density of dust forming el- 
ements residing in some not closer specified dust components, 
and (2) since it is assumed that growth of dust is essentially re- 
stricted to molecular clouds (cf. Draine |T"9901 l which relates the 
dust growth problem, like the dust destruction problem, closely 
to the multiphase structure of the ISM, which has approximately 
to be taken into account. 

It is generally assumed that the growth of dust grains of a 
specific kind / is governed by some rate determining reaction 
step, usually by the addition of that one of the elements required 
to form the chemical compound which has the lowest abundance 
in the gas phase, and that the rate of addition of all other more 
abundant elements adapt to the slowest process. The growth is 
determined in this case by some specific key element and a spe- 
cial atomic or molecular species from the gas phase carrying 
most or all of this key element, the growth species. The key ele- 
ments for the condensed phases of interest are given in Table [5] 
The equation for the change of the mass mj of a single grain of 
species j is 
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(25) 



Here n Jgr is the particle density of the growth species, v ; , t h, g r its 
thermal velocity, a ; the growth coefficient, J[ the surface area 
of the grain and Aj the atomic weight of one formula unit of 
the dust material under consideration. v ; m and v/, c are the num- d f \ 
ber of atoms of the key element contained in the growth species 
and in the formula unit of the condensed phase, respectively. 
Evaporation is neglected since this is not important at the low 
temperatures in molecular clouds. The change of the mass den- 
sity pj of the dust particles in a molecular cloud is obtained by 
multiplying the growth equation of single grains by the number 
density of grains and the probability distribution of grain radii 
(assuming spherical grains) and integrating over all grain radii a. 
One obtains 



Fig. 12. Growth timescale for the dust species growing in molec- 
ular clouds for the Milky Way model at the solar cycle: Silicate 
dust (full line), carbon dust (dashed line), iron dust (dotted line) 



key element, possibly lowered by the fraction of this element 
that is blocked in some unreactive molecular species. Let / de- 
note the fraction of the key element already bound in the dust 
species under consideration, the degree of condensation f. The 
gas-phase density of the growth species is (1 - /) n^max- Hence 
we obtain the following equation for the degree of condensation 
in the molecular cloud 
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Numerically we have 
46 Myr x 
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where V\j = Ajm AMV /p c is the volume of one formula unit in the 
condensed phase, p c is the mass density of the condensed phase, 
and (. . .) denotes averaging with respect to the size distribution 
of grains. We define the following average grain radius 
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(a) 

For a MRN size distribution (Mathis, Rumpl & Nordsiek [T977l l 
we have for instance 



(a), - V fl oa l ~ 0.035 pm , 



(28) 



where ciq = 0.005 pm and a\ = 0.25 pm are the lower and upper 
limits of the distribution of grain radii, respectively. 

The maximum possible particle density of the growth species 

is 



where T,; gr is evaluated with characteristic values for the physi- 
cal variables. The temperature of clouds is assumed to be 10 K, 
the growth coefficient a at such low temperatures is assumed 
to be a — 1. The characteristic growth time is generally short 
compared to the lifetime of molecular clouds, except at very low 
metallicity of the ISM. As an example Fig.[T2lshows the growth 
timescale t, gr calculated from our Milky Way model at the solar 
cycle for the important dust species. 

In principle, the average grain radius (a) 3 depends on the 
degree of condensation / (a oc f 1 ^ for compact structures), but 
we neglect this weak dependence. In this case the equation for / 
can immediately be integrated with the result 
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(32) 



where A^h is the number density of H nuclei in the molecular 
cloud (usually equal to 2n^), e the element abundance of the 



Here /o is the initial degree of condensation at t = 0. For t » Tgr 
the degree of condensation approaches / = 1 . 
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4.2.2. Source term for dust production 

Molecular clouds form in the interstellar medium by instabili- 
ties, mainly during the compression of ISM material in the snow- 
plow phase of SN shocks. They disappear within a rather short 
time if active star formation starts and winds of massive stars 
and expanding supernova bubbles disperse the clouds. For the 
average lifetime of molecular clouds we take an observationally 
and theoretically motivated value of « 1 X 10 7 yrs (Leisawitz 
et al. [19891 Will iams & Mc Kee [19971 Matzner |2U02l Krumholz 
& McKee 120061 Blitz et al. 120071 This val ue for t he lifetime is 
somewhat shorter than that used in Tielens ( 1998) in his model 
of dust growth in clouds, but seems to be more appropriate for 
the most massive clouds, which contain nearly all of the ISM 
mass in clouds. The lifetimes of the clouds equals the charac- 
teristic timescale T exc h by which matter is exchanged between 
clouds and the remaining ISM. 

At the instant of cloud formation the clouds inherit the dust 
content of the interstellar medium outside of clouds. The dust 
content of the matter outside of dense clouds is lower than within 
clouds since dust destruction processes operate in this material, 
while in clouds the dust grows by accreting not yet condensed 
refractory elements. In fact, except if the metallicity of the ISM 
is very low, the growth timescale is much shorter than the life- 
time of the cloud and the condensation of the refractory elements 
runs into completion before the cloud disappears. 

Let the initial degree of condensation of the key element for 
some dust species be /o. If after a period t a cloud is rapidly 
dispersed, the degree of condensation in the matter returned to 
the ISM material outside clouds is equal to the value given by 
Eq. ( T32l >. The effective dust mass return for species j by a molec- 
ular cloud is then 
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(34) 



(1 +4e He )v/,c 

is the maximum possible mass-fraction of the dust species in the 
material of the molecular cloud and M c ioud is the cloud mass. In 
principle one has to observe that some fraction of the cloud mass 
is converted into stars and not converted into other phases of the 
ISM. Since we describe the effect of dust consumption by star 
formation within the frame of our approximation by a separate 
term in Eq. (0, this process needs not to be accounted for in 
Eq. Q3). 

Equation ([331 has to be multiplied by the probability P(f) that 
the cloud is destroyed at some instant within the period between 
t and t + dt 

P(t) = _L e -' /r »* , (35) 

^"exch 

and integrated over t. Here it is assumed that the cloud destruc- 
tion occurs at random with a mean lifetime r exc h. Finally, averag- 
ing with respect to the mass spectrum of clouds and multiplying 
with the surface number density of clouds, one obtains for the 
mass return rate of MC-grown dust per unit time and unit area 
of the galactic disk 

Gj,d = (//.ret _ /j,o) -^/.max^cloud > (36) 
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where £ c i ou d is the surface mass density of clouds and the average 
degree of condensation on cloud dispersal is 
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(37) 



The quantity G^j is the gain term that has to be used in Eq. (0 
for the evolution of the MC-grown dust component j. 

In principle the evaluation of this term requires to consider a 
multiphase ISM where molecular clouds form one of the com- 
ponents. Since we wish to consider the simpler model of a one- 
phase ISM we have to cast Eq. (l36l l to an appropriate form for 
this case. In terms of the mass fraction of clouds in the ISM 
Xcioud = S c i oud /S ISM we have 
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is the maximum possible surface density of dust of kind j if 
all material from the ISM that can be condensed into this dust 
species is really condensed, and Xy jd is the mass fraction of dust 
of kind j in that part of the ISM that is not in clouds. 2, is the 
surface mass density of the key element for species j in the ISM. 
For the average mass-fraction of dust in the total ISM we have 
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which yields 
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It follows 
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We define the effective exchange time 
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This is much longer than T exc h since X c ioud 1 an d reflects the 
fact that it requires many timescales T exc h to cycle all ISM ma- 
terial through clouds where it is laden with fresh dust. Our final 
result for the dust production term is 
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(42) 



This is the appropriate dust production term for MC-grown dust 
in the approximation of a one-phase ISM model. In this model 
the mass-fraction X c i ou ^ of the ISM in clouds is a free parameter 
which does not follow from the model calculation but has to be 
taken from observations. We will use a value of X c i ou d = 0.2 
(cf. Tielens 2005), which is appropriate for the ISM at the solar 
cycle. 

4.2.3. Limit cases 

The degree of condensation / ; ret in the material returned from 
clouds at the time of their dispersal depends essentially on the 
ratio of the growth timescale r gr to the average cloud lifetime 
Texch- If Tgr s> r exc h (slow growth at low metallicities) one ex- 
pects that only small amounts of dust are added to the initial dust 
content, in the opposite case (rapid growth at normal metallici- 
ties) one expects complete condensation in the returned material. 
This can be confirmed by calculating the lowest order terms of a 
series expansion of the integral in Eq. (f37j for the two limiting 
cases. 
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T exch ' T gr 



Fig. 13. Approximation for the variation of the degree of con- 
densation / let with Tgr/Texch for fo = 0.3. The full line shows the 
result of a numerical evaluation of the integral (l37l i. the dashed 
lines the two limit cases Eqs. d43l and d44b . and the dotted line 
the approximation I 



For slow growth (j sr » T exc h) one introduces f/T exC h as inte- 
gration variable, expands exp [ (T exc h/T gr ) t ] in a series, and inte- 
grates term-by-term. The result is in the linear approximation 
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(43) 



If this is inserted into Eq. ( 1331 one recognizes that T exc h cancels 
out and, thus, the amount of dust produced during the residence 
time of ISM material in the cloud phase is determined by the 
details of the growth process. 

For rapid growth (r exc h s> T gr ) one introduces f/r gr as inte- 
gration variable, expands exp [ (T gr /T exc h) t ] in a series, and inte- 
grates term-by-term. The result is in the linear approximation 
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(44) 



(note that ln(l - fo) < 0). If this is inserted into Eq. ( f33l > one 
finds that the dust production by the clouds is nearly indepen- 
dent of the details of the growth process within the clouds and 
is (almost) exclusively determined by the cycling frequency of 
ISM material between the clouds and the other phases of the 
ISM. The composition of the dust, of cause, is also in this case 
determined by the details of growth processes. 

The variation of f J)iet with T gl -/T excri in both limit cases for a 
value of fo = 0.3 are shown in Fig. [13] as dashed lines together 
with the result of a numerical evaluation of the integral (|37| >. that 
is shown as full line. A rather accurate analytic fit formula for 
the full range of T gr /T excrl values is 
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This approximation is shown in Fig.Q~3]as dotted line. The re- 
sults for other values of fo are similar. Only for very small fo the 
approximation becomes somewhat worse in the transition region 
Tgr/Texch ~ 1, but for bigger fo it becomes even better. For the 
purpose of model calculations it suffices to use the approxima- 
tion d45i 



4.3. The individual dust species 

Evaluation of the source term Eq. ( f42l > for dust requires to cal- 
culate the growth time scale t, gr , given by Eq. (|30l , Xj >rmK given 
by Eq. (l34b . and the degree of condensation /}, ret in the returned 
material, which we calculate from the approximation d45l ). for all 
dust species j which are formed by growth in molecular clouds. 

The constants required for calculating these quantities are 
given in Table [5] The growth coefficient is assumed to be a = 1 
for all cases since at the low temperatures in dense molecular 
clouds of about 10 K even the weak attractive van der Waals 
forces lead to adsorption; the basic theory for this is discussed, 
e.g., in Hollenbach & Salpeter (119701 1, and Watson (|1975l l. 

For calculating the average (0)3 we use in all cases the 
approximation Eq. (f28b following from a MRN-size distribu- 
tion (Mathis et al. 119771 1. This is only a crude approximation, 
but without attempting to calculate grain size distributions it is 
hardly possible to fix this quantity with better accuracy. 

The initial value fj$ for calculating f^ ret is given by the de- 
gree of condensation in that part of the ISM matter which is not 
in clouds, i.e., one has 
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Using this and Eq. d45T > in Eq. ( f39b yields in principle a non- 
linear equation that has to be solved for fjo- For most pur- 
poses it suffices to replace Eq. d46b by the approximation / 7> o ~ 



Xj,d/Xj tB 
big. 



since the difference between Xj& and Xm is not very 



4.3.1. Silicates 

The silicate dust in the ISM accounts for about one half of the 
total dust mass (e.g. Dwek 2005), but its composition is still a 
matter of debate. Studies of silicate composition based on in- 
terstellar depletions, modelling of extinction curve, and in situ 
measurements of dust in the local ISM give quite different re- 
sults, although they all agree on olivine ([Mg x Fei_ x ]2SiC>4 with 
< x < 1) and pyroxene (Mg x Fe]_ x Si03 with < x < 1) 
as major candidates for ISM silicates. A number of studies of 
depletions of Mg, Fe, Si and O atoms in the interstellar gas 
phase came to the conclusion that observed depletions indicate 
an olivine-type stoichiometry of dust in the diffuse ISM (Savage 
& Sembach 1996, Jones 2000). In a recent attempt to fit the sili- 
cate features of the interstellar extinction curve Min et al. (2007) 
found that the composition of the ISM silicates is consistent with 
a Mg-rich mixture of olivine and pyroxene with a bigger contri- 
bution from pyroxene than from olivine. Fitting of the 9 and 18 
fim features of the extinction curve shows that, while the 9 /mi 
feature can be fitted well by olivine dust, the position and peak 
strength of 18 /mi feature is fitted much better with a pyroxene- 
type stoichiometry (Demvk [1999l l. 

An olivine-pyroxene mixture with a contribution of more py- 
roxene than olivine is therefore chosen for modelling the ISM 
silicates in the present paper. As a first approximation we adopt 
a fixed silicate composition to study silicate dust production by 
dust growth in molecular clouds. Modelling of a variable sil- 
icate composition, depending on local growth conditions, is a 
challenging problem to be considered in future papers. Let f Q \ 
be the (fixed) fraction of the silicate dust that has olivine stoi- 
chiometry; the fraction 1 - f \ then has pyroxene stoichiometry. 
Assuming the same Mg fraction x for both olivine and pyroxene 
in our model, two parameters determine the silicate properties: 
f Q \ and x. 
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The total efficiency of dust production by molecular clouds 
does not show a significant dependence on the choice of the 
parameters f \ and x. Variations of the Mg-fraction x change 
the total dust mass on the level of 10% at most, but define the 
silicate-to-iron dust mass ratio. This is due to the fact, that for 
the Mg-rich mixtures that are considered here, Mg is the critical 
growth species. With decreasing x, less Mg is needed for sili- 
cate dust growth, but the total silicate mass increases due to a 
bigger contribution from the Fe-bearing component while at the 
same time less Fe remains for growth of solid iron. We fix the 
Mg fraction x to a value of x = 0.8 by fitting the present-day 
silicate-to-carbon dust mass ratio of the model to its observed 
value of 0.6, inferred from observations of the infrared emission 
from the Diffuse Infrared Background Experiment (Dwek et al. 
[19971 . 

The olivine fraction f D \ is chosen such as to reproduce the 
observed Mg/Si ratio in dust using the simple relation for a given 
olivine-pyroxene mixture: 

Ajvig 

xA si 



fol = 
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(47) 



Here AM g and As; are observed abundances for the elements Si 
and Mg bound in dust (in particles per million hydrogen atoms, 
ppm). The ratio AM g /Asi equals 1.06 or 1.07 for dust in the 
diffuse ISM, as given by Dwek d2005l or Whittet ([20031 . re- 
spectively, which results in a value of f \ = 0.32. Although the 
Ajyig/Asi ratio obviously varies in different ISM phases, we use 
average dust abundances from the diffuse medium, since this 
constitutes a significant fraction of the total ISM mass, and only 
little is known about the very cold dust in molecular clouds. 
Test calculations for different AM g /Asi ratios available from dif- 
fuse ISM studies showed no strong influence on dust masses, 
resulting in 4% change of total dust mass with 10% decrease of 
A M g/Asi ratio. 

For given silicate composition, the growth species used 
to calculate the growth timescale, Eq. ( [301 , is determined by 
the abundance of the least abundant species available for dust 
growth. This is either Si or Mg and we chose in Eq. d30l 
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where vsi < 



v Mg,c 



1.06. 



4.3.2. Carbon dust 

The formula unit is the C atom, i.e., one has v, c = 1. It is as- 
sumed that C is present in the gas phase in molecular clouds pre- 
dominantly as free atoms or in a number of molecules bearing 
one C atom only and that these serve as growth species. Some 
fraction fco of the carbon is blocked in the CO molecule and is 
not available for carbon growth. The precise fraction cannot be 
fixed without calculating models for the chemistry of the molec- 
ular clouds. Observations indicate a CO abundance in molecular 
clouds of 20% . . . 40% of the C abundance (e.g. Irvine et al. 
[19871 van Dishoek et al. [19931 van Dishoek & Blake [19981 . 
In the calculation we consider the two cases £co = 0.2 and 
£co = 0.4. The carbon abundance e in Eqs. (|30l and (|34l is 
calculated as 

£c - ^C.sic 
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Fig. 14. Growth of dust in molecular clouds at the solar cycle. 
Thick lines show /, the average degree of condensation of the 
key elements into dust for the dust species shown at the instant 
when the molecular clouds are dispersed and their material is 
mixed with the other phases of the ISM. Thin lines show /o, the 
corresponding degree of condensation at the formation time of 
clouds. One always has /o < / since dust grains grow in molec- 
ular clouds and are partially destroyed again in the ISM out- 
side of clouds until they enter the next cloud. Growth of iron 
dust in clouds starts with a significant time delay because of de- 
layed iron production by SN la events. The calculation is for 
£co = 0.2; the result for £co = 0.4 is not shown because the 
corresponding curves are almost the sames. 



where Sc, s ic is the surface density of C bound in silicon carbide 
dust 



v _ 12 



4.3.3. Iron dust 
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1. It is as- 



The formula unit is the Fe atom, i.e., one has 
sumed that Fe is present in the gas phase as free atoms, which are 
the growth species in this case. The iron abundance e in Eqs. ( [30l 
and ( 1341 is calculated as 
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56S H 

where 2 Fe s j] is the surface density of Fe bound in silicate dust 
species. 

5. Results of the evolution model for the solar cycle 

5. 1 . Evolution of MC-grown dust 

The model for the dust evolution considers silicate dust, car- 
bon dust, and iron dust as species that grow in dense molecular 
clouds. The corresponding growth timescales r gr calculated from 
our Milky Way model at the solar cycle are shown in Fig. [12] 

During the first Gyr of evolution of the galactic disk the 
metallicity at the solar cycle is low ([Fe/H]^ -2, cf. Fig. [3J and 
the characteristic growth timescale of dust in clouds exceeds the 
average lifetime of dense molecular clouds of about 10 Myr as- 
sumed in our model. Only small amounts of dust are added to the 
dust content of the interstellar matter during its cycling through 
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Fig. 15. Evolution of the dust mass fraction in the in- 
terstellar medium of the main interstellar dust compo- 
nents and of the Stardust species at the solar cycle. The 
dust grown in molecular clouds dominates the total dust 
mass of the interstellar medium. For carbon dust two 
results are shown corresponding to an assumed frac- 
tion of 0.2 resp. 0.4 of the carbon in molecular clouds 
blocked in the CO molecule. The dust condensed in 
stellar ejecta (AGB stars, SNe) has an only small abun- 
dance in the ISM. The condensation efficiencies of dust 
in supemovae used for the model calculation are given 
in Tab. [5] 



clouds. This can be seen in Fig.[14]that shows the evolution of 
the initial value fj$ of the degree of condensation of the key ele- 
ments into dust, defined by Eq. d46b . for each of the dust species 
j, and the average final degrees of condensation fj, calculated ac- 
cording to Eq. (f45b for the same species, if the clouds are finally 
dissolved. Both quantities, fjo and /;, are calculated during the 
course of our model calculation for the evolution of the Milky 
Way at the solar cycle. During the first about one Gyr the degree 
of condensation of refractory elements into dust increases only 
marginally by dust growth in molecular clouds. Therefore, the 
dust production in the Milky Way is almost completely deter- 
mined by dust condensation in the ejecta of stars, and the dust 
content of the ISM is determined during this transient phase by 
dust injection from stars into the ISM and by dust destruction in 
the warm phase of the interstellar medium. Obviously the devel- 
opment would be considerably different if one has a strong star- 
burst at early times and metallicity becomes already high before 
the first AGB stars appear, but this seems not to have happened 
in the case of our Milky Way. 

Once the metallicity of the ISM has grown to a level of about 
[Fe/H]= -2 some dust starts to condense during the lifetime of 



molecular clouds and their dust content at the instant of their dis- 
solution somewhat exceeds their initial dust content. From this 
point on molecular clouds start to contribute to dust production 
in the galaxy. 

If the metallicity has climbed after more than 2 Gyrs to a 
level of about [Fe/H]= -1, the degrees of condensation into dust 
fj at cloud dispersal are much higher than the degrees of conden- 
sation into dust fjo at cloud formation; in fact, dust growth runs 
almost into completion during the lifetime of the clouds. During 
each cycling step of interstellar matter through clouds the matter 
is laden with fresh dust and this dust is mixed into the general 
ISM at cloud dispersal. The dust content of the ISM then is es- 
sentially determined by the equilibrium between dust growth in 
clouds and dust destruction in the warm phase of the interstellar 
medium. 

The degree of condensation / of carbon into carbon dust 
does not approach unity (see Fig. [T4l . since it is assumed that 
20 to 40% of the carbon in molecular clouds forms CO and then 
is no more available for dust condensation. 

The iron dust abundance evolves somewhat different from 
that of the silicate and carbon dust. The main reason is that most 
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of the Fe is produced in SN la explosions and these turn on rather 
late due to the long lifetime of their low mass precursor stars. 
Additionally we assumed in our model that SN la explosions do 
not start until the metallicity of the precursor stars has rised to 
[Fe/H] k, — 1 (see Sect. 12.1.5b . A second reason is that in our 
model for dust growth it is assumed that the silicates grown in 
clouds contain a certain fraction of iron and the small fraction of 
iron initially produced by supernovae is then almost completely 
consumed by the growth of silicates with some iron content. This 
will somewhat change if the iron content of the silicates is not 
fixed, as in our present calculation, but will be determined from 
growth kinetics. 

5.2. Evolution of dust abundances 

Figure Q3] shows the evolution of the various dust components 
during the 13Gyrs of evolution of the galactic disk. The dust 
components with index TSM' are the isotopically normal grains 
grown in the interstellar medium. Surviving grains from stellar 
sources are characterized by an index AGB' or 'SN' if they are 
from AGB-stars or from supernova ejecta, respectively. 

The results depend on the efficiency of dust production by 
stars, dust condensation in molecular clouds, and dust destruc- 
tion rates in the interstellar medium. The dust production by low 
and intermediate mass stars on the AGB is determined from the 
table of Ferrarotti & Gail (2006 ) and the dust destruction rate 
from Jones et al. (1996). They are probably not too far from 
reality. The dust production efficiencies of massive stars are un- 
known. One has, however, one piece of information: the abun- 
dance ratios of the presolar dust grains from AGB stars and SNe. 
We have varied the supernova dust production efficiencies 77 in 
Eqs. ( TT6b . . . ( fT9l until the observed abundance ratios for silicate, 
carbon, and SiC dust from AGB and SNe sources is reproduced. 
Details are described in Sect. 15.31 the resulting efficiencies are 
listed in Table[5] These efficiencies are very low, probably since 
they account also for a number of destruction effects that prevent 
dust formed in SNe from escaping into the general ISM. 

The dust population of the ISM in this model is dominated 
by dust grown in molecular clouds except for the very earliest 
times, where Stardust dominates. The model shows that preso- 
lar dust grains with their isotopic anomalies revealing the origin 
of these grains are always a minor component of the interstellar 
dust. Most of the dust in the ISM has collected nearly all of its 
material from the interstellar gas phase and is isotopically incon- 
spicuous. If new stars are formed from the ISM containing such 
a dust mixture, the dust in their protoplanetary accretion disks 
contains an only tiny fraction of presolar dust grains with iso- 
topic anomalies. This fits well with the recent findings obtained 
with the nano-SIMS investigations of Interplanetary Dust Grains 
(IDPs) by Messenger et al. (|2p03 ), which show that nearly all of 
the silicate grains from cometary nuclei, which should be domi- 
nated by interstellar grains, are isotopically norma0 

The population of Stardust grains is dominated by grains 
from AGB stars because of the low efficiency of SN dust pro- 
duction. In our model the AGB dust is dominated by carbon dust; 
silicate dust and SiC dust is much less abundant. In meteorites 
presolar carbon dust in the state of graphite is the least abundant 
of these three components (cf. Nguyen et al. [2007 )■ The dis- 
crepancy is certainly due (i) to the different survival properties 
of different kinds of dust material in the Solar System and the 



5 It is a little bit puzzling that the STARDUST particles analyzed so 
far seem to be mainly material from the Solar System (see Zolensky et 
al. |2006l McKeegan et al. [2006b 



Table 6. SN presolar dust fractions and corresponding derived 
efficiencies of dust production. 



Dust species 


Fraction of SN grains 
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0.01 


0.00015 


Silicates 


0.03 


0.00035 




0.05 


0.0005 




0.1 


0.04 


Carbon 


0.3 


0.15 




0.5 


0.20 


SiC 


0.01 


0.0003 



parent bodies of the meteorites, and (ii) the methods of labora- 
tory investigations applied for different dust grains. This disables 
presently a comparison between abundances of different presolar 
species predicted by the model and observed in meteorites. 

One outstanding feature of the abundance evolution of preso- 
lar dust grains is the rather late appearance of silicate and SiC as 
compared to carbon grains. This reflects the fact that AGB stars 
synthesize the carbon required for soot formation from He and 
have not to rely on external sources of heavy elements. Contrary 
to this, the Si bearing dust components cannot be formed until 
sufficient Si is synthesized in supernova explosions and returned 
to the ISM from which subsequent stellar generations inherit the 
Si required for formation of Si-bearing species. This needs some 
time and additionally the precursor stars of the main sources of 
Si-bearing dust, the AGB stars, are rather long lived low-mass 
stars (cf. Fig. [Toll. Presolar silicate dust grains in the ISM where 
a rather new phenomenon at the instant of Solar System forma- 
tion. 

The low abundance of silicate Stardust may also provide an 
explanation for the lack of crystalline silicate dust in the ISM 
though a lot of crystalline dust is injected to the interstellar 
medium by outflows from AGB stars. Even if there would be no 
amorphization processes with energetic electrons and ions (cf. 
Demyk et al. 2004; Jager et al. 120031 1. crystalline silicate dust 
(< 20% of the silicate dust injected by AGB stars) would be too 
rare compared to amorphous ISM dust to be observable by its 
absorption features. 

5.3. Efficiency of supernova dust production 

The dust mass produced in the ejecta of supernovae is not 
known. Observations indicate that only small amounts of dust 
condense and that only part of all SNe form dust. With the kind 
of model for dust evolution in the ISM we have developed, one 
can try to determine for some dust species an estimate of the ef- 
ficiency of dust production by supernovae. This can be done by 
comparing the abundance ratios of supernova dust and AGB dust 
resulting from the model calculation with real observed abun- 
dance ratios of presolar dust grains with SN and AGB origin 
in meteorites. Only silicate, carbon, and SiC dust are presently 
suited for this, because only for these dust species the required 
data for presolar dust grains are available. Iron dust has not yet 
been detected as presolar dust so far and it is unclear whether it 
really exists. 

Such a comparison depends on some assumptions. The first 
one is that the production rate of dust by AGB stars is known 
with significantly better accuracy than the dust production rate 
of supernovae. The second basic assumption of a comparison be- 
tween these kinds of data is that the fraction of the dust destroyed 
between the instant of its incorporation into the just forming 
Solar System and the instant of laboratory investigation of preso- 
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lar dust grains does not depend on the kind of stellar sources 
where the dust has formed, but only on its chemical composition. 
One has to assume, in other words, that the basic properties of 
AGB and SN dust with the same composition, its size spectrum, 
its resilience against oxidation and treatment by strong acids, are 
the same. 



5.3.1. Silicon carbide 

The observed abundance ratio of X-type SiC grains and 'main- 
stream' SiC grains in the presolar dust population isolated from 
meteorites is close to 0.01 (Hoppe et al. 2000). Fitting the ef- 
ficiency 77si C ,SNii in Eq. < TT~8T > such that the calculated abundance 
ratio for SiC from supernovae of type II and AGB stars agrees 
with the observed abundance ratio yields ?7sic,SNii = 5 x 10~ 4 . 
This ratio seems surprisingly low, but already the low abundance 
of X-type SiC grains shows that the efficiency of SiC dust for- 
mation in supernovae is low. The efficiency 77 s i c ,SNii of SiC dust 
condensation in supernova determined in this way is used for the 
final model calculation and is that given in Table[5] 

The abundance ratio for the SiC grains refers to grain abun- 
dances observed after isolating the grains from the meteorite ma- 
trix by a rather brutal treatment with oxidizing agencies and by 
strong acids (cf. Amari et al. I1994K but it has been argued by 
Amari et al. (119941 1 19951 1 that at most a small fraction of the 
grain material is lost during this procedure. On the other hand, 
the size distribution of SiC grains in Murchison meteorite found 
by Daulton et al. (2003 ) shows a lack of grains smaller than 0.5 
yum diameter that dominate in circumstellar dust shells (e.g. Jura 
1 19971 1, i.e., the small grains are lost already in the ISM or in the 
Solar System. If there should be severe systematic differences in 
the mass fraction of sub-micron sized grains in the size distribu- 
tions of SiC grains of SN and AGB origin, the abundance ratio 
derived from isolated SiC grains may be severely misleading, 
but presently we have no better data. 



5.3.2. Silicate dust 

The number of silicate grains from stellar sources detected in 
meteorites and interplanetary dust particles is small up to now 
(Nguyen et al. 2007; Messenger et al. 2005 ). Besides the about 
some 100 silicate grains with isotopic anomalies attributable to 
an origin from AGB stars only a single grain has been detected 
up to now with isotopic characteristics pointing unambigously to 
a SN origin (Messenger et al. 2005 ). A few more have been de- 
tected that are likely also of SN origin (cf. Vollmer et al. 120071 ). 
The small numbers do not allow to pin down the abundance ratio 
of silicate dust from the two possible sources with any reliabil- 
ity, but as a working hypothesis we assume an abundance ratio 
of 3%. Then we can determine an efficiency of silicate dust for- 
mation in supernovae of 77 s ii.sNii = 3.5 x 10~ 4 . This is the value 
given in Table [5] The efficiencies for a somewhat smaller and 
bigger abundance ratio are shown for comparison in Table [6] 

Contrary to the case of SiC the silicate grains are detected 
by scanning technics from material that has not been prepared by 
chemical treatment. There is therefore no problem to be expected 
in the sense that part of the grain population is already destroyed 
by preparation methods before the particles are investigated. 

5.3.3. Carbon dust 

The abundance of presolar graphites from supernovae is highly 
uncertain. Chemically separated graphite fractions were further 
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Fig. 16. Evolution of abundances in dust per million of hydrogen 
atoms of the main dust-forming elements as predicted by the 
model calculation. For carbon two lines are shown. The upper 
one is for the case that a fraction of £co = 0.2 of the carbon is 
bound in the inreactive molecule CO, the lower one for the case 
£co = 0.4. 



subdivided into low density separates KE1 and high-density sep- 
arates KFA1, KFB1 and KFC1 (Amari et al. [l994l While many 
- though not all - high density graphites seem to have an AGB 
star origin (Croat et al., 2005 ), low density graphites are ascribed 
to supernovae (Hoppe et al. 1995 Amari et al. 1995; Travaglio 
et al. 1999 ), particularly inferred from isotope data of the low 
density fraction KE3 (Amari et al. |1995bl ) which is the coarse 
grained (> 2 jjm) subgroup making up 70% of KE1. 

If all low density graphites are from supernovae, this 
would correspond to a relative abundance of 67% (by weight). 
However, there are significant uncertainties which fractions of 
the various density separates do indeed correspond to a specific 
supernova or AGB star signature, so we adopt here an abundance 
of 50 ± 30 % (Hoppe, pers. comm.) and calculate 3 different 
cases for 10%, 30% and 50% of all graphites coming from su- 
pernovae. For the model results shown in the figures we assumed 
a mass-fraction of 30%. From this one derives an efficiency of 
carbon dust formation in supernovae of 77 s ii,SN n = 0.15. This is 
the value given in Table[5] This efficiency is much higher than in 
the two preceding cases and would mean that in SN are mainly 
sources of carbon dust. Efficiencies for a smaller (10%) and big- 
ger (50%) abundance ratio are shown in Table|6]for comparison. 

Presolar graphite grains have mainly size ^ 1 yum (e.g. Zinner 
119971 ) while for carbon dust grains in circumstellar dust shells 
around AGB stars one knows that they have sizes 0.1 fim. 
Only the small fraction of grains from a large-size-tail of the size 
distribution are found in the separates investigated in the labora- 
tory. If the graphite grains formed in SN ejecta should have sys- 
tematically bigger sizes than that formed in AGB-star outflows 
(there is however no indication for this), the supernova graphite 
dust fraction found in the separates would overestimate the true 
abundance of graphite grains from supernovae and our estimated 
efficiency t/s^snii would be too high. 

On the other hand, a much higher condensation efficiency 
of carbon dust seems not unlikely since then only the carbon 
atoms in the carbon layer have to condense into dust particles 
and no complicated mixing processes of the supernova ejecta 
are required as for the other dust species. 
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Fig. 17. Predicted average depletions of main dust forming ele- 
ments at present time are shown with filled circles. Upper and 
lower points for depletions for carbon calculated with CO mass 
fraction 0.4 and 0.2 correspondingly. Upper and lower error bars 
represent observed depletions in warm and cold diffuse clouds 
from Welty et al. ( fl999l for C, Si, Fe and Cartledge et al. d2006l 
for O and Mg. Filled triangles marks the average depletions in 
diffuse clouds (see Whittet 2003 and references therein). 



5.3.4. Need for dust accretion in the ISM 

The low efficiency of dust production by supernovae indicated 
by the rather low abundance of Stardust of SN origin compared 
to Stardust from AGB stars means that the supernovae cannot 
contribute substantially to dust in the ISM, contrary to what is 
frequently assumed. Therefore it is unavoidable that most of the 
dust mass observed to exist in the ISM is formed in the ISM 
itself and not in stars. This has consequences for the dust pro- 
duction in young galaxies with low metallicity, where only su- 
pernovae can be sources for Stardust. The high dust abundances 
observed in some high-redshift galaxies can, according to our 
results, not result from the first generation of SNe but requires 
already additional accretion processes of heavy elements in in- 
terstellar clouds. This kind of dust evolution in young starburst 
galaxies will be treated in a separate paper. 

5.4. Depletion of elements 

Fig- US] shows our model results for the abundance evolution of 
the main dust forming elements. Since we assumed a fixed com- 
position of ISM silicates, the ratio between Mg, Si, and O does 
not change during evolution, but this is not the case for Fe, which 
is consumed both by silicate and iron dust production. In the fig- 
ure two lines are shown for carbon, corresponding to two differ- 
ent assumed fractions £co of the carbon in molecular clouds tied 
up in the inreactive CO molecule. The upper one corresponds to 
£co = 0.2, the lower one to £co = 0.4, bracketing typical ob- 
served values (e.g. van Dishoek & Blake [T998j l. A higher value 
of £co means that less carbon is available for dust formation. 

The observed gas phase abundances of elements in diffuse 
interstellar clouds indicates various degrees of depletions of 
many of the dust forming elements relative to their solar abun- 
dances. This is explained as resulting from their condensation 
in interstellar dust. The amount of the dust forming elements 
locked up in interstellar dust (this is what is shown in Fig. [Tol l, 
however, cannot reliably be derived from observations. The stan- 
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Fig. 18. Evolution of the dust-to-gas ratio at the solar cycle as 
predicted by the model calculation. 



dard procedure is to determine instead gas phase abundances of 
the elements and subtract these from some kind of 'standard' 
cosmic element abundances in order to determine how much of 
each element is condensed into interstellar dust (cf. Sembach & 
Savage 1996 for a review). However, to draw conclusions about 
the dust composition from observed depletion patterns one needs 
to make a decision of what set of abundances is used as the ref- 
erence abundances for the elements. Frequently Solar System 
abundances, or abundances of nearby F & G stars or of B stars 
are adopted (cf. Tables [2] and |3), resulting in different dust com- 
positions. 

A modelling of the chemical evolution of the Galaxy includ- 
ing dust allows a study of the evolution of the depletion of the 
gas abundances by dust condensation and a comparison of the 
model with present observed data, since gas and dust abundances 
are known from calculations. However, a one phase ISM model 
reflects properties of the dust averaged over the different ISM 
phases, therefore only a qualitative comparison of depletions 
is possible. One should notice that observed depletions are re- 
stricted to diffuse clouds, while molecular clouds are too opaque 
to be studied in absorption lines. 

Our predicted averaged depletions, at present time and at the 
solar cycle in the ISM, for the 5 main dust forming elements un- 
der consideration (C, O, Mg, Si, and Fe) are shown in Fig. [T7] 
For comparison, observed depletions in warm and cold diffuse 
clouds from Welty et al. d 19991 1 and Cartledge et al. (2006), and 
average depletions in diffuse clouds (see Whittet 2003) also are 
shown in the figure. The model calculation reasonably repro- 
duces the observed values, except for a somewhat low calculated 
degree of depletion of Fe. 

The degree of depletion of iron cannot be increased by as- 
signing a much longer destruction timescale t ; snr for iron dust. 
A model calculation shows that this does not significantly in- 
crease the depletion because the lifetime of dust grains is in 
any case limited by the timescale of dust consumption by star 
formation, which is about 2.3 Gyr and hence already not really 
big compared to the lifetime against destruction by shocks. Also 
a higher than assumed stability of Fe-bearing silicate does not 
help, since then a higher than observed depletion of Si is to be 
expected. The main reason for the low degree of depletion in 
the model seems to be that the fraction of Fe in the gas phase 
is not completely determined by the destruction of Fe-bearing 
grains by shocks in the warm phase but to a significant extent 
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Fig. 19. Composition of the interstellar mixture of dust species grown in molecular clouds, the MC-grown dust, and of the presolar 
dust species from AGB stars and supernovae, the Stardust, at the solar cycle. Left: At the instant of Solar System Formation. Right: 
For the present solar neighbourhood 



also by return of gas-phase Fe by stars, which needs some time 
until it is cycled into clouds and depleted from the gas phase 
by dust growth processes. In our model the degree of depletion 
of Fe (and this holds in principle for all refractory elements) is 
limited by the rather long time required for cycling of matter 
between the ISM matter not in dense clouds and the matter in 
dense clouds. The origin of the low degree of depletion of Fe 
in the model calculation is presently unclear, but probably bears 
physical significance and may indicate that for iron some slow 
accretion of Fe atoms from the gas phase into dust is possible 
also in the warm and/or cold phase of the ISM. 

5.5. Evolution of the dust to gas ratio 

Figure [18] shows the evolution of the dust-to-gas ratio according 
to the model calculation. The hydrogen gas-to-dust mass ratio is 
approximately 100 for the diffuse ISM averaged over long lines 
of sight passing through a number of interstellar clouds (Spitzer 
119541 . Recent studies of the hydrogen gas-to-dust ratio in the 
Local Interstellar Cloud (Kimura et al. 2003a) also confirms the 
canonical value from Spitzer d!9541 l. If one converts this to a 
ratio of dust mass to total gas mass, one gets a value of 0.007. 
The value of the dust-to-gas ratio in our model for the present 
time ISM is close to the average value derived from observations. 
Our model therefore nicely reproduces the average dust mass 
fraction of the Milky Way in the solar vicinity. 

5.6. Dust input into the Solar System 

Figure [19] shows the composition of the interstellar dust mix- 
ture at the instant of Solar System formation and the present day 
composition. Numerical values are given in Table [7] Both mix- 
tures are not significantly different since the abundances of re- 
fractory elements in the ISM have changed only slightly over the 
last 4.56 Gyr (cf. Fig. [5]). This dust mixture is clearly dominated 
by MC-grown dust and contains an only small fraction of star- 
dust. The Stardust is dominated by dust grains from AGB stars; 
dust grains with SN origin form only a minor component. The 
composition seems to be consistent with models for the compo- 



Table 7. Surface densities of different dust species at instant of 
Solar System formation and at present time as predicted by the 
model calculation. 



Dust species 


Source 


£;,d( f SSF)[M pC 2 ] 


%d(f G )[M oP c- 2 ] 




AGB 


5.39-10- 3 


1.1210 4 


Silicates 


SN II 


2.51-10" 6 


2.08- 10~ 6 




ISM 


2.76- 1(T 2 


2.76- 10~ 2 




AGB 


7.73-KT 4 


5.73-10~ 4 


Carbon 


SN II 


7.86- 1(T 5 


7.73- 10~ 5 




ISM 


1.43 10 2 


1.5410 2 




AGB 


1.8910 3 


2.75-10- 5 


Iron 


SN II 
SN la 


7.72- 10- 7 
9.83- 10- 7 


6.63- 10- 7 
9.55-10- 7 




ISM 


5.31-10 3 


6.78-10- 3 


SiC 


AGB 


3.80- lO" 5 


4.58- 10" 5 


SN II 


8.2M0~ 7 


6.82-10- 7 



sition of the local interstellar dust (Kimura et al. I2003bl Frisch 
l2006l l 

The dust mixture at time of Solar System formation is that 
from which the solid bodies in our planetary system formed. 
Relics of this dust mixture can be found in the Solar System in 
two types of objects: matrix material of primitive meteorites and 
in comets. However, until Solar System bodies form from the 
dust component of the matter collapsed from some part of the 
parent molecular cloud into the protoplanetary accretion disk, 
the material underwent a number of alteration processes. Even 
the most primitive material in Solar System bodies is not just 
unmodified ISM matter. For this reason meteoritic matrix ma- 
terial is presently not suited for a comparison with the model 
results, since the alteration processes on the parent bodies are 
presently not completely understood (cf. McSween et al. 2002). 
Material from comets is probably more suited, and once more 
detailed results from the STARDUST mission are available, it 
may be possible to compare the model predictions for the ISM 
dust composition entering the protoplanetary accretion disk of 
the Solar System with observations. Presently most of the ana- 
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lyzed particles returned by the STARDUST mission seem to be 
material from the solar system (Zolensky et al.[2006; McKeegan 
etal.HXj6]l. 

Presently one can only state that the dust mixture inherited 
by the Solar System from its parent molecular cloud and the 
present day ISM dust mixture in the solar neighbourhood pre- 
dicted by our model calculation are roughly in accord with the 
dust composition estimated by Pollack et al. (|1994| l from obser- 
vations of extinction properties of the dust material and consid- 
erations on element abundances, which is presently held for the 
best estimate of the composition of the dust material from which 
the Solar System formed. 



6. Concluding remarks 

In the present paper we have developed a model for calculating 
the chemical evolution of of the Milky Way disk and the dust 
content of the interstellar medium in a consistent fashion, fol- 
lowing partially the methods proposed by Tielens ( 11998b and in 
particular by Dwek ( I1998I ). 

The chemical evolution part of the model for the Milky 
Ways disk follows mostly well-established methods. The model 
is checked with the standard tests applied to such type of models 
and reasonably reproduces the observational constraints. A new 
aspect is the use of the new tables of Nomoto et al. (2006]! for 
the heavy element production by massive stars. This gives better 
agreement between the calculated evolution of the element abun- 
dances of the main dust forming elements (O, Mg, Si, Fe) and the 
observed evolution of abundances in the ISM as witnessed by the 
variation of atmospheric element abundances of main sequence 
G stars with metallicity. In particular the problem with the low 
Mg abundances disappears with the new results of Nomoto et 
al. (120061 I. A good reproduction of the abundance variations of 
the dust forming elements is important if one tries to model the 
interstellar dust mixture. 

The model for the evolution of the interstellar dust is based 
on three essential elements: 

First, the model uses for the dust input to the ISM by low and 
intermediate mass stars the results of the model calculations of 
Ferrarotti & Gail (2006), which combine synthetic AGB evolu- 
tion models with models for circumstellar dust shells including 
dust formation in the stellar wind. These models describe for 
the first time consistently the dependence of dust production by 
AGB stars on stellar initial mass and metallicity. 

Second, the dust production by supernovae is described by 
a simple parametrization, that was already applied by Dwek 
(1998). A lot of information has accumulated over the years 
by the efforts of the meteoritic science community on studying 
nucleosynthetic processes in stars by analyzing isotopic abun- 
dances in presolar dust grains. As a by-product this provides us 
with abundance ratios of presolar dust grains from supernovae 
and from AGB stars. This allows for the first time to estimate 
the efficiency of dust production by supernovae by fitting cal- 
culated abundance ratios of Stardust in the Solar System to ob- 
served abundance ratios of presolar dust grains from AGB stars 
and SNe. Applying these gauged efficiencies yields the unex- 
pected result, that dust production by massive stars is not impor- 
tant for the evolution of the ISM dust component, at least not for 
Pop II and Pop I metallicities. 

Third, a simple approach is developed to include dust growth 
in molecular clouds into a model for the evolution of the inter- 
stellar dust. This approach is in a sense incomplete, since the 
growth of dust (as its destruction) depend heavily on the phase 



structure of the interstellar gas (hot, warm, cold), which can- 
not be treated adequately within the approach generally used 
in chemical evolution models. So some quantities in the dust 
growth model like the mass fraction of ISM matter in clouds 
have to be taken from observations, but we have data for these 
quantities only for the present day Milky Way. Since the distribu- 
tion of the ISM over the phases is neither spatially nor temporal 
constant, in a realistic modelling the phase structure should be 
part of the model calculation. Nevertheless the present approach 
allows for the first time to study the evolution of the interstellar 
dust population, including several Stardust species. 

The results obtained with this model are in reasonable accord 
with, for instance, the observations of interstellar depletions of 
refractory elements. Again, one has the problem that the degree 
of depletion depends on the phases of the ISM, which are not 
adequately considered in the present model. Presently we have 
started work to implement our dust model in a chemodynamical 
evolution code (Berczik et al. 120031 ) which removes the short- 
comings of the present model. 
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